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FOREWORD 


This report describes a nonlinear rotorcraft me del and associated 
computer software which has been developed and documented for NASA, Langley 
Research Center, Hampton, Virginia under contract NASl-1^570 (July 1976). 
This work has been performed by the Lockheed-California Company, Burbank, 
California. 

P. H. Kretsinger (Lockheed) performed the software implementation. 

W. D. Anderson and Fox Conner (both of Lockheed) assisted in preparation of 
the program. 
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REXOR II ROTORCRAFT SIMULATION MODEL* 

Volume I — Engineering Documentation 
J. S. Reaser and P. H. Kretsinger 
Lockheed-Califbrnia Company 
SUMMARY 

This report describes a generalized format rocorcraft nonlinear 
simulation called REXOR II. The program models single main rotor vehicles 
with up to seven main rote" blades. Wings, two horizontal tail planes, 
and auxiliary thrustors may be included to model a variety of compound 
helicopter configurations. 

Program output is primarily in the form of machine plotted time his- 
tories specified from a signal list. This list is, in turn, user selected 
from a set of computation variables used by the program. 

LIST OF SYMBOLS 


The symbols used in the REXOR II equations are quite numerous. In order to 
keep the catalog of symbols to manageable proportions the following list 
is divided according to the discussion in Section 3. Namely, a list of 
basic symbols is given, followed by subscripts, superscripts, and post- 
scripts. Nonconforming cases of usage together with complicated or 
obscure subscripting are fully annotated in the basic list. 

SYMBOLS 


u. 

a 

s 



1 

[A] 

*1,2,3 



arbitrary vector 
speed of sound, m/s 

2 

acceleration vector, m/s (ft/s ) 

longitudinal component of blade first harmonic flapping, vad. 
generalized mass element matrix 
modal variables 

generalized displacement of nth blade, first mode 


•The contract research effort which has lead to the results in this report 
was financially supported by USARTL (AVRADCOM) Structures Laboratory. 



2n 

A 3n 

*18 

b 

B 

B 1S 

c 

Cd 


"M 


X,Y,Z 

C l,2,3 

c(k) 


d 

dr 

dt 

d/dt 

(d/e) ( 

(d/e)., 

e 

El 

f iMF 


generalized displacement of nth blade, second node 
generalized displacement of nth blade, third node 
cosine component of blade first harmonic cyclic, rad 
number of main rotor blades; arbitrary vector 
dissipation function 

sine component of blade first harmonic cyclic, rad 

blade segment chord, m (ft) 

damping matrix • 

aerodynamic drag coefficient 

aerodynamic lift coefficient 

aerodynamic pitching moment coefficient 

power coefficient 

thrust coefficient 

linear damping, N/m/s (lb/ft/s) 

rotary damping, N-m/rad/s (ft-lb /rad/s) 

blade bending to feathering couplings 

lift deficiency function 

infinitesimal increment 

increment in rotor, radius, m (ft) 

increment in time 

derivative with respect to time 

swashplate to feather gear ratio, zero collective 

swashplate to feather gear ratio slope with collective 

pitch horn effective crank arm, m (ft) 

blade bending stiffness distribution, N-m 2 (lb- ft 2 ) 

ground effect factor for main rotor 
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F 

F 

X.Y.Z 

F 

'M.'J' 

F 

BPH 

e 

g X,Y,Z 

G 

{ 0 } 

G 

GJ 


XX 


YY 


ZZ 


XY 


XZ 


YZ 


J 

J 

k 


W 


factor; force, N (lb) 

force components along X,Y,Z directions, N (lb) 
generalized force about <fi, 6, axis 

feathering mode generalized force 

2 2 
gravity, m/s (ft/s ) 

gravity components along X,Y,Z directions 
gear ratio 

generalized force vector 

gyro angular acceleration partial product 

2 2 

blade torsional stiffness, N-m (lb - ft ) 

= Enu X^ 2 , kg-m 2 (slug- ft 2 ) 

= Enu Yj 2 , kg-m 2 (slug- ft 2 ) 

= Enu Z^ 2 , kg-m 2 (slug- ft 2 ) 

= Em i (Y ± 2 + Z i 2 ), kg-m 2 (slug-ft 2 ) 

= Em i (X i 2 ♦ Z i 2 ), kg-m 2 (slug-ft 2 ) 

= Em i (X 2 + Y i 2 ), kg-m 2 (slug-ft 2 ) 

= Enu X^ Y^, kg-m 2 (r; lug- ft 2 ) 

= Em^ X^ Z^, kg-m 2 (slug-ft 2 ) 

= Em^ Y^ Z^, kg-m 2 (slug-ft 2 ) 
unit vector 
unit vector 
advance ratio 

number of blade radial stations; reduced frequency, 
rad/s; unit vector 

spring matrix 

blade spring matrix element 


3 



4b 


■f 

■h 

■i 

^S» 

N 

tM] 

"rk 

"x 

"y 

*Z 

*X.Y,Z 

*♦ 

I 

P 

P iMR 


spring constants aloog X,Y,Z direction, 1/a (lb /ft) 

spring rates about ♦, 6, ♦ axis, S-a/rad (ft-lb /rad) 

location inboard feather bearing, a (ft) 

location outboard feather bearing, a (ft) 

radial Iccatioo of intersection of precoae and feather 
axis, a (ft) 

rolling aoaent, K-a (ft- lb) 

aass of eleaent, kg (slugs) 

sinned fuselage coordinate aass, kg (slugs) 

saned hub axis aass, kg (slugs) 

aass of i th particle or blade se^co*, kg (slugs) 

svashplate sinned aass, kg (slugs) 

pitching aoaeat, M-m ( ft- lb ) ; « Ea^, kg (slugs); aach nuaber 
generalized aass eatrix 
generalized aass Matrix eleaent 

* E *i X i* **"■ ^ slu 8 -ft ) 

* laj Y it kg-a (slug-ft) 

■ Ea i Z^ kg-a (slug-ft) 

a oa en ts about X,Y,Z axis, 1-a (ft-lb) 
blade torsional aoaent, H-a/a (ft-lb /ft) 
nuaber of system particles 

angular velocity about X axis, rad/s; particle 

main rotor pitch aoaent inflow, m/s (ft/s) 

generalized coordinate; angular velocity about Y axis, 
rad/s 
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4 iMR 

Q 

q a 

QLOADS 





t 

T 


CT3 

u 

U 

U %P,S,T 

v 


w.,_ 

iMR 

v 

iTR 

x 


main rotor roll moment inflow, m/s (ft/s) 
generalized forcing function 

aerodynamic pressure times reference wing area, kg (lb) 
total nonmain rotor aerodynamic loads matrix 
tail rotor torque, N-m ( ft- lb ) 

general vector; radius of curvature, ft; angular velocity 
about r axis, rad/sec; notation for (X,Y,Z) 

static blade shape 

vector displacement of particle p in X,Y,Z axis system 

vector displacement of x,y,z origin in X,Y,Z system 

gyro damper coupling ratios 

Laplace variable, path of motion of particle p 

blade spline length along neutral axis locii, a (ft) 

time 

kinetic energy, H-m ( ft- lb ) 
transformation of coordinates matrix 
velocity in X direction, m/s (ft/s) 

potential energy function, N-m (ft-lb); strain energy, N-m (ft-lb) 

air velocity on blade element , m/s (ft/sec) 

velocity in Y direction, m/s (ft/sec) 

trajectory velocity 

velocity in Z direction, m/s (ft/sec) 

main rotor collective inflow, m/s (ft /sec) 

tail rotor collective inflow, m/s (ft/sec) 

motion in X direction, m (ft); blade span location 
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X 


y 

x 


TTO 


1,2.3 


OHA 


z 

Z 


SP 

Z TTO 


1,2,3 


Z 


OBL 


OF 


OSP 


$ 


0 


FA 

B PHn 


coordinate direction; axis; deflection, a (ft); location, 
a (ft); cross product 

blade radial station of sweep and Jog, a (ft) 
trajectory path, a (ft) 
tail rotor longitudinal force, a (lb) 
notion in ¥ direction, a (ft) 

coordinate direction; axis; deflection, a (ft); location, a (ft) 
tension torsion pack, outboard end nodal coefficients 

difference between Y direction locations of eg and neutral 
axis points of blade eleaent, a (ft) 

notion in Z direction 

coordinate direction; axis; deflection, a (ft); location, a (ft) 

relative swash? late vertical displacement with respect to 
the hub, a (ft) 

tension- torsion pack outboard end nodal coefficients 

teetering rotor undersling, a (ft) 

hub set distance above fuselage set, a (ft) 

hub set distance above svashplate set, a (ft) 

angle of attack, rad 

angle of attack with hub set, rad 

sideslip angle, rad 

blade feathering angle, rad 

feathering/pitch-horn bending or dynanic torsion 
generalized coordinate displaceaent 

blade droop relative to precone angle, zed 
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y blade sweep angle* rad; dynamic stall delay, s 

Y ^ trajectory path angle with E set, rad 

6 limit deflection, rad; freeplay, rad; s m all increment 

5 tail rotor pitch - flap coupling 

jTn 

dc/da downwash factor of wing on horizontal tail 

C vector notation of +, 8, 

6 rotation about Y axis, rad 

0q collective blade angle, rad 

A sideslip at blade element, rad 

p air density, kg/m\ (slugs/ft^) 


T 


T 


0 


♦ 



X iMR 

* 


* 


c 


^PH 


* 


T 


* 


W 


u 


time constant, s; natural period, s 
feathering axis precone, rad 
rotation about X axis, rad 
feathering angle, rad 

feathering angle of blade element of nth blade, rad 

blade roc:, reference feather angle, rad 

blade torsion, rad 

sum of blade twist and torsion, rad 

wake angle of main rotor, rad, (deg) 

rotation about Z axis, rad; sideslip angle with hub set, rad 

control input axis rotation from swashplate, rad 

pitch lead angle, rad, (deg) 

trajectory path yaw with E set, rad 

main rotor apparent airflow angle, rad 

rotational speed, rad/s; angular velocity, rad/s; 
natural frequency, rad/s 
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partial derivative, derivation 


3 

SUBSCRIPTS 

a 

A 

b 

HERD 

BLE 

BLn 

C 

CG 

CORR 

DU 

DIM 

E 

EHG 

EST 

F 

FA 

FB 

Fn 

FR 

G 

GEN 

GFB 

G3P 


arbitrary coordinate set a 

due to aerodynamics 

arbitrary coordinate set b 

associated with blade elastic bending 

blade element coordinate system 

blade reference axis system for the nth blade 

associated with pilot control input, chordvise 

associated vith center of gravity location 

corrective, correction 

referring to dovnvash 

referring to dynamic component 

earth axis 

associated vith poverplant - engine 
estimated 

fuselage axis; associated vith blade feathering 
referring to blade feather axis 
associated vith feedback 

associated vith feathering of the n th blade 
due to friction 

referring to gyro or gyro coordinate system 
associated vith gas generator section of poverplant 
associated vith gyro control feedback 
gyro to svashplate connection 
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GUB relating to gyro gimbal unbalance 

K referring to hub or principal reference axis system 

HT associated with horizontal tail 

i referring to inflow, particle 

IB referring to inboard feather bearing location 

j spring matrix index 

jog associated with blade attachment Joggle 

J associated with gyro end of feedback rod linkage 


Jn 

k 

LAG 

LIMIT 

m 

MR 

n 

NA 

HEW 

NO 

NR 

OB 

OLD 

P 

PH 

r 

R 


associated with feedback rod coming from the n th blade 

generalized mass index 

associated with lead-lag damper 

signifying limiting value 

blade mode index, spring matrix index 

associated with main rotor 

blade number index 

referring to blade segment neutral axis 

newly determined value 

normal (to airflow) component 

pertaining to nonrotating value 

referring to outboard feather bearing location 

value from previous time step 

associated with propeller; perpendicular blade component 
referring to pitch horn 
generalized mass index 
referring to rotor axis system 
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REF 

RM 

S 

SC 

SP 

SP 

c 

S, SP 
STEADY 

SV 

T 

TR 

TRIM 

TV 

UB 

UNSTEADY 

VT 

WING 

X 

Y 

YA 

Z 

ZA 


associated with blade feather reference value 

referring to control gyro feedback lever moment 

referring to blade spanvise velocity; general mode; static; 
structural; shaft 

referring to blade segment shear center 

referring to swashplate 

comnand to svasbplate 

referring to svashplate limit stop 

steady component 

referring to blade sweep angle location 

associated with trajectory path relating to E axis; 
tangential blade component; blade torsion; blade twist 

associated with the tail rotor 

initial or trim value 

associated with blade twist (built in) 

relating to control gyro unbalance 

associated with unsteady component 

associated with vertical tail 

associated with the wing 

relating to component in X direction 

relating to component in Y direction 

relating to aerodynamic component in Y direction 

relating to component in Z direction 

relating to aerodynamic component in Y direction 
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0 


1,2,3 

IS 

1/4 c 
3/4 c 

S PHn 

* 

6 

* 

SUPERSCRIPTS 

I 

T 

(-) 

(’) 

(•) 

(") 

(- 1 ) 

(-) 

POSTSCRIPTS 

(i) 

(n) 


(nought) associated with collective value, coordinate axis 
value, with respect to principal reference axis, blade 
root summation 

with respect to blade modes 1, 2, or 3 
first harmonic component shaft axis feathering 
with respect to blade 1/4 chord 
with respect to blade 3/4 chord 

associated with the feathering mode of the n th blade 
relating to component in the <J> direction 
relating to component in the 0 direction 
relating to component in the direction 

referring to inertial reference 

matrix transpose 

(bar) average quantity 

(prime) slope with respect to blade span 

(dot) time derivative of basic quantity 

(double dot) second time derivative 

matrix inverse 

vector quantity 

blade radial station index 
blade number index 
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1. INTRODUCTION 


1.1 Scope of the REXOR II Program 

REXOR II is a rotorcraft analysis tool which has resulted from applying an 
interdisciplinary math modeling philosophy. The REXOR II math model is 
written for a single rotor helicopter with capability for analysis of hinge- 
less or hinged rotor systems vith conventional controls. This helicopter 
may be conventional in design, winged, or compounded. The main rotor may 
have a maximum of seven blades. The model is broken dowr. into the three 
major categories shown in Figure 1. These categories are the control system, 
the rotor, and the body. 

Figure 1 indicates the manner in which these components are related to 
one another as utilized in the f ’-'•sis. The analysis is the simulation of 
an entire aircraft, which indue detailed dynamic description of the 
rotor and control system as well . conventional six-degree-of-freedom 
body dynamic description which operates in two modes identified as TRIM 
and FLY. In the TRIM mode, the aircraft is constrained to a prescribed 
static flight condition while the controls are activated and the rotor is 
allowed to respond to obtain a force and moment equilibrium of the aircraft 
at that static condition. In the FLY r.ode the entire aircraft is free to 
respond dynamically to control inputs or to any other arbitrary inputs such 
as gusts. Pilot inputs can be any single or multiple control manipulation 
in the form of simple steps or pulses, doublets, stick stirs, or other 
transient input within the capabilities of the control system simulated. 

As a result, transient loads and resulting aircraft and rotor dynamic re- 
sponse can be obtained. For correlation purposes, actual flight test con- 
trol motions can be used as input to provide comparative response data. 
Additionally, gust inputs and other types of external excitations can be 
applied directly to the rotor and/or airframe. 


1.2 REXOR II Capabilities 

REXOR II is a detailed rotorcraft math model simulation with particular 
emphasis on the main rotor mechanics. The program is particularly valu- 
able in a detailed exploration of rotor characteristics of proposed de- 
signs, in identifying problem areas and verifying fixes in flight test 
development programs. A case history is given in Reference 1. 
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Figure 1. - Block diagram model description 



Typical REXOR II applications are listed below. 


Dynamics : 

• Rotor stability at. a function of flight speed, maneuvers, 
rotor rpm, nonlinear blade aerodynamics 

• Rotor/body sensitivity and dissipation capacity as a function 
of gusts and oilot control inputs 

• Effects of design parameters (mechanical and elastic 
couplings, controls, etc.) on rotor stability and load 
sensitivity 

• Correlation and check of specialized dynamic models. 

Handling Qualities : 

• Vehicle response to pilot control inputs for vehicle flight 
conditions, speed, altitude, rotor rpm, design parameter 
variations 

• Vehicle stability as function of speed, rotor rpm, flight 
conditions, design parameters 

• Effect of design parameter variations on handling qualities 

• Development and checking of handling qualities models. 

Failure Analysis : 

• Effect of loss of one inplane damper on subsequent flight 
time history 

• Blade projectile hit and ensuing events 

• Blade strike and resulting rotc-r track. 

Performance : 

• Correlation and independent check of performance models, 
particularly in regions of highly nonlinear blade aerodynamic 
operation (retreating blade stall and compressibility effects) 

• Develop data for performance models for use in nonlinear 
areas 


ll+ 



Loads : 


Steady-state rotor loads as a function of rotor rpm, flight 
velocity, control trim settings 


• Dynamic rotor loads as a function of rotor rpm, flight velo- 
city, vehicle maneuvers, pilot control inputs 

• Rotor /fuselage clearances as a function of speed, vehicle 
maneuvers, rotor rpm, pilot control inputs, flight 
configuration 

• Rotor/ fuselage /wing design characteristics requirements as 
functions of maneuver load factor, control commands (see 
Reference 2). 


1.3 Improper Application REXOR II 

While REXOR II is capable of performing a number of analysis tasks, the 
progrrm range of use is certainly not all inclusive. Examples of types 
of use where REXOR II either wouldn't work well or would be impractical are 
given below. 

REXOR II is an extensive math model and, as such, may consume a considerable 
amount of computer time to execute a case. Therefore, the program is not 
intended as a parametric design analysis tool, but rather as a device to 
verify the correctness of a parametric selection process. 

REXOR II does not treat blade-to-blade vortex interaction. This condition 
limits the validity of the vibration solution in the transition flight 
regime. 

REXOR II typically uses tventy or less blade radial stations. The computer 
blade deflections show good correlation to measured data with this model- 
ing. However, since shear is a first derivative, and moment is a second 
derivative of deflection data, care needs to be exercised in their use 
(Reference 3). 


1.1* The REXOR II Report and Its Use 
This report is presented in three volumes. 

• Volume I 

A development of rotorcraft mechanics and aerodynamics including 
a derivation of the equations of motion from first principles. 
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Volume II 


The development and explanation of the computer code required to 
implement the equations of motion. 

• Volume III 

A user's manual containing a description of code input/output end 
instructions to operate the program. 

Volume I is intended to be a self-suf 'icient guide to the math development 
of the equations of motion and is the reference background as such. Volume 
II gives the location of computation elements , and serves to locate elements 
for inspection or modification. Volume III presents normal program operation 
plus troubleshooting gu'de material required for day-to-day program use. 
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2 


BASIC COMPUTATIONAL IDEA 


2.1 Modal Solution - Overview 

The aircraft is described dynamically by an array of fully-coupled degrees of 
freedom. In addition to the six degrees of freedom of the fuselage principal 
reference axes, six degrees of freedom describe rotor hub to fuselage deflection 
due to shaft bending and transmiss5.cn mount motion. Rotor/engine speed is a 
degree of freedom. The control svashplate has three degrees of freedom. 

Motion of each of the main rotor blades is described by three coupled flap- 
wise and inplane modes and a pitch horn bending degree of freedom which couples 
blade feathering to the swashplate. The total number of degrees of freedom 
possible is 16 + ub, where b is the number of blades. 

The blade modes are primitive modes in that they are determined from a 
lumped parameter analysis at a select: d rotor speed and collective blade 
angle, hereafter referred to as the reference feather angle. The general- 
ized stiffness matrix is computed using these rotating modes and contains 
only the structural stiffness of the blades and hub. This formulation 
ensures proper internal and external force and moment balance. The model 
deflections outboard of the feather hinge are rotated through the actual 
feather angle less the reference feather angle. Thus, blade element 
deflections outboard of the feathering hinge due to modal displacements are 
defined to remain aligned with a coordinate axis system which is orthogonal 
to a plane containing the instaneous deformed feather axis and rotated through 
the instantaneous feather angle less the reference f .at her angle. As a result, 
the internal strain energy in the blade due to unit model displacements is in- 
variant with variation in blade angle. This technique permits the highest 
resolution of motion and forces for the blade with an assumed mode solution 
for a given number of modes. 


2.2 Ehergy Methods Development 

The equations of motion for REXOR II are developed from Lagrange's equations, 
which is an energy approach. If one can express the kinetic, potential, and 
dissipative energies of a system in addition to the work done by external 
^orces, then Lagrange's equations provide a powerful method for developing 
he equations of motion. 

The dynamic equations of motion are written in matrix form as 


- [A] |q| + | Gf = 0 


(1) 
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where [ A ] is a square matrix of generalized mass elements, )q } is a column 
vector of accelerations of the generalized coordinates and J G f is a column 
vector derived from the Lagrangian energy functions, dissipation function and 
generalized forces, which take the form: 




The equations cf motion are solved using a time history solution with rotor 
azimuth angle increments required to provide a stable solution at the highest 
frequency mode present. 


2.3 Calculation of Rotor Mode Displacements, Velocities, and Accelerations 

In a rotor simulation of this type, it is difficult to compute the proper dis- 
placement velocities and accelerations and associated inertia and aerodynamic 
forces and moments which are required for high resolution of the blade feather- 
ing moments. This requires exacting aerodynamic data as well as a precise 
statement of the inertial loadings. To establish the feathering moments due 
to these loads, the relationship between the x'eather axis and the point of ap- 
plication of the loads must be precisely determined. This is accomplished by 
a very accurate analytic construction of the undeformed blade and a superposi- 
tion of the blade elastic bending on this shape. In order to achieve the 
highest resolution of the predicted blade shape and feather axis position, the 
blade modes are defined at approximately the trim collective blade angle. The 
blade static position is also constructed at this blade angle. Blade element 
displacements, velocities, and accelerations are then computed from the com- 
bined static shape, the elestic blade motion, and blade feathering with re- 
spect to the reference feather angle. 

The aerodynamic description used in the analysis is composed of a rotor inflow 
model, nonlinear steady and unsteady blade element aerodynamics, nonlinear 
fuselage aerodynamic characteristics, rotor /body aerodynamic interference, and 
auxiliary airloads from the tail rotor and tail surfaces. The main rotor down- 
wash effect on the wing and horizontal tail angles of attack is an empirical 
function of rotor thrust and forward velocity. The nonlinear fuselage aerody- 
namics nay be inputted as tables of actual wind tunnel test data. 

The aircraft primary control systems are simulated from the pilot control 
levers operating through a boost system in all control axes. Gearing and 
gains in the control path are inputs to the analyst:, and nriy be easily changed 
for studying the effects of design changes in the control system. 

Control servos are simulated by first-order lags with rate limits and with 
soft and hard physical stops. Contrl stiffnesses in collective and cyclic 
pitch axes of the main rotor are included in the dynamic equations of motion. 
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2 U Output 


The analysis is a time history solution of the equations of motion. REXOR II 
does not directly process the results of the solution process for output , it 
creates an output file of user selected parameters which are correlated by 
the computation time step. Prom this data bank the recorded signals can be 
selected for tabular or plotted output. Assuming a good selection of param- 
eters is chosen to be recorded, the user in an interactive mode may select 
as little or as much of the information for viewing as is needed. Thus a 
configuration can be examined thoroughly without having to rerun the case 
to select additional output. 
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3. SYMBOLS 


The notation used in REXOH II generally follows what could he tensed RASA 

notation. In general: 

• Axis systems use a right-hand triad X, Y, Z 

• Rotations about these axes are also a right-hand triad 8 » *, $ 

• Rotation rates, again a right-hand triad,are p, q, r 

• Velocity components of X, Y, Z are u, v, w. 

3.1 Subscripting Rotation 

Subscripting is used as a rule in REXOR II to further identify a variable. 
Superscripts except in a few column vectors are reserved to denote raising 
to a power. The subscripting can mean: 

• IVpe of element; F for fuselage, SP fcr swashplate, TR for 
tail rotor, R for rotor, etc. 

9 Coordinate system reference; BLn for blade axis, H for hub axis, 
R for rotor axis, etc. 

• Modal identifiers. 

3.1.1 Blade number . - The blade medal identifier typically is of the form 

A . Where n is the blade number, 
am 

3.1.2 Mode number- - Also from A y m is the mode number, and is keyed to 

1 ■ mn 

the symbol A. A represents blade bending modes (3). Therefore m can be 
1 to 3. 

3-1.3 Mode type . - Other than blade bending the remaining blade mode is 
torsion, and is separately identified as . Nonblade nodes are identi- 
fied by the direction and subscripted axis of notion. Examples are for 

n 

rotation of the rotor and e g for shaft pitching. 
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3 -l- 1 * Generalized mass, dapper, spring, forces . - The generalized masses 
are denoted as M doubly subscripted by the two modes active for that mass. 
Examples are M and M. _ . This scheme is also used for other elements 
♦S*S A mn ®H 

of the equations of motion, dampers (C), springs (K), forces (F). Note the 
forces are a column vector and singly subscripted. 

3-1.5 Forces and moments . - In the process of forming the equations of 
motion many subelements of forces and moments are formed, translated and 
combined. Several levels of subscripting may exist in performing this 
process. The guidelines to the layering are: 

• First level denotes the direction or axis system that the quantity 
is formed in. Examples are X and BLE. 

• Second is the axis system involved or axis system being translated 
to, depending on the specification of the first level. The second 
level may also be specified as 0 or nought, to indicate the 
value is a', the coordinate system origin. This notation is used 
to show an inertial reference and blade root summation quantities. 

• The third level, usually outside a series of bracketed quantities , 
shows the blade number being computed, or the overall coordinate 
system in use for the computation at hand. 
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U. COORDINATE SYSTEMS AND TRANSFORMATIONS 


1* . 1 Introduction 

Prior to developing the equations of motion, a system of coordinate sets 
with a description of the elements of the system in these sets and the 
interrelationship of the sets is required. 


h.2 Coordinate Sets 

1.2.1 Fuselage coordinates {Xp, Yp, Zp). ” The fuselage X and Z axes lie in 
the fuselage plane of symmetry. The location of the origin is arbitrary. See 
Figure 2. The coordinates form a right-hand triad Xp, Yp, Zp. Notations for 
velocities with respect to earth of these coordinates are either Xp, Y p , Zp 
or Up, Vp, Vp. A conventional double dot notation is used for acceleration. 
Euler rotations of the set follow conventional practice of roll right $p, 
pitch up 8p, and yaw right ♦p. Rates of rotation are either denoted by dot 
notation or Pp, qp, r^,. Angular acceleration is double dot notation of the 
rotation or dot notation of the rates, ip, 6p, or Pp, 4p, ry. 

Numerous aerodynamic terms are referenced to the fuselage set. Figure 3 
shows the relationship of airflow to this set. The components of airflow, 
also noted as Up, Vp, Wp, are defined with respect to the fuselage set by 

an angle of attack a, and a sideslip angle 8. The angle of attack is the 
arcsin of the ratio of the vertical component and the vector sum of the X 
and Z components. The sideslip is the Y component of airflow in relation 
to the total vector airflow sum. The angle of attack is positive (pitch up) 
of the fuselage set with respect to the airflow. The sideslip is positive 
(yaw left) for the airflow relative to the set. The airflow is the vector 
sum of the fuselage set inertial motion and flow fields from other parts of 
the vehicle, such as main rotor downwash. 

1.2.2 Hub coordinates (Xg 

the main rotor mast, but does net rotate with the mast. 

Airflow information is referenced to the hub set for use in the main rotor 
aerodynamic calculations. The reference scheme is shown on Figure 1. For 
components of airflow u , v , w with respect to the hub set, an angle of 

4i n n 


Yg, Zg). - The hub set origin is at the top of 
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a = 



-1 


\/» f 2 + * f 2 4 


Figure 3. - Coordinate systems fuselage axis 





attack a^, and sideslip \p are defined. The generatici conventions are different 
from the* 1 fuselage airflow reference in order to clearly separate the inplane 
and outplane airflow components. 

-.2.3 Shaft coordinates (X,, Y_, Z_). - The shaft axes are an intermediate 
■ s s s 

set between the hub and fuselage sets, see Figure 5- The geometry is deter- 
mined by the distances (X , Y , Z )_ „ and (X , Y ,Z )_ „ that the hub origin 

O o or-n o oo b-n 

is located from the fuselage and shaft axes and by the rotation (<J> , 6 , 

o o 

^ )„ of the shaft axes from the fuselage axes. The hub axes are parallel 

o S— K 

to the shaft axes. 

The elastic deflections due to motions of the shaft and transmission suspen- 
sion are given by the set of coordinates (X,Y Z , * , $ , y, ) . The hub 

s s s s s s 

is assumed to move as a rigid body with respect to the shaft axis origin. 

U.2.U Rotor coordinates (Xp, Yp, Zp). - The undeflected rotor set has the 
same origin as the hub set. See Figure 6. The Xp and Yp axes rotate with 
the blade number 1 reference axis system. At 4<p = 0 

aligned but point in a direction opposite to the X^ and Z axes. The rota- 
tion of the rotor set is measured counterclockwise (CCV) from the -X^ axis 
by the angle ♦ p. 

L.2.5 Blade coordinates (X_ T , Y__ , Z„ ) . - To hookkeep the deflections 
BLn BT.n BLn 

properly for all the main rotor blades, sets equivalent co the rotor set 
are created for each blade. These are the BLn sets, where n is the blade 
number (counted clockwise from blade number one). All BLn sets are iden- 
tical except for an azimuthal rotation (n - 1) A where is the inter- 

blade angular spacing. The rotation is about the Z axis. Note that BLn 

R 

sets are rotating coordinates and have a common Z axis. 

L.2.6 Blade element coordinates (X„ TT ,, Y^ T ^, Z„ TT ,). - The blade element set 
" BEE BLE BLE 

origin is located at the center of gravity of an element of a particular 

blade. See Figure 7. Reference to a column vector subscripted by BLE is 

used to denote the blade element located by the blade element set origin. 

The right-hand coordinate triad of this set has the X axis parallel to the 

local quarter chord line, the Y axis along the chord line toward the leading 

edge. The Z axis is mutually perpendicular and pointed up. The BLE set is 

used to track the local feather angle, to develop aerodynamic and dynamic 

loading terms. 


the 


and 




axes are 
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Figure 5- - Coordinate systems - hub, shaft, and fuselage sets. 
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Figure 7. - Coordinate systems - blade element set. 


The BLE set origin for each blade element specifies the element c.g. with 
respect to the quarter chord, and in terms of the BLE directions, i.e., 
for the Kth element the position coordinates are SY(K) and SX(K) where 
SX(K) is the blade radial station. Transformations to the neutral, no- 
stretch axis are made for X deflections. Note: The quarter chord is 

merely a convenient reference datum, and does not convey any model 
limitations or assumptions. 

I*. 2. 7 Freestream (earth) set (X„, Y^, . - The freestream set is essen- 

E E E 

tially the earth or inertial set inasmuch as the axis alignments are the 
seme. However, the freestream set can assume any origin. Thus the use of 
the set is to reference the local gravity vector and/or an absolute angular 
displacement or linear velocity acceleration of another set. As shown on 
Figure 8 , the Zg axis points down toward local gravity. Other sets refer- 
ence to the E set, as the F set shown here, may assume any starting value 
of roll and pitch such as the trim initial conditions. The relative ori- 
entation changes with progressing time of flight. 


With the freestream set origin located coincident with the fuselage set, 
the components of fuselage set velocity in E set are Ug, Vg, Wg. These 
components combine into a trajectory velocity u^ and path X^,. The 
trajectory path is yawed right and pitched up y t from the E set. 

See Figure 9. 


U.2.8 Swashplate coordinates (Xg p , Y gp , Z ). - As shown on Figure 10, the 

SP set origin is located in line with the Z axis and above the hub set a 

ri 


distance Z 


OSP' 


The SP set does not rotate with the rotor shaft . For no 


deflection of the SP set, the X and Y axes have the same alignment as the 
X and Y of the hub set. 


b.3 Degrees of Freedom 

The degrees of freedom of the REXOR II equations are defined as the general- 
ized coordinate variables of the set of equations of motion to be devel- 
oped in Section 5. These degrees of freedom fully describe the motion of 
*he physical elements of the modeled helicopter, but each direction of 
otion of the helicopter may not have a degree: of freedom directly asso- 
ciated with it. The physical motions may be described by a series of 
modal variables (Section 5.U) or through a set of transformations and com- 
binations of the degrees of freedom as developed in Sections b . U ar.d U . 5 . 
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Figure 8. - Coordinate systems - freestream (earth) 
to principal reference axis. 






t H (RIGHT) 


Figure 10. - Svashplate coordinate system- 



The REXOR II rotorcraft simulation analysis can be applied to describe the 
vehicle-rotor -control system dynamic response for up to 16 + l*b (where b is 
the number of blades) fully-coupled degrees of freedom. These include the 
normal six rigid body or vehicle degrees of freedom; rotor speed; provisions 
for up to twenty-eight degrees of freedom defining rotor blade motion (four 
mode degrees of freedom for seven blade maximum) . Three swashplate degrees- 
of-freedom and six for describing shaft /transmission deflection. The equa- 
tions of motion are written in a general form so that additional degrees-of- 
freedora can be added if desired. The current degrees -of -freedom are listed 
in Figure 11. The discussion following describes them in detail. 


U.3.1 Vehicle or rigid body . - The six rigid body degrees-of-freedora; three 
translations, and three rotations, are defined as motions of the fuselage or 
principal reference axis system. Section U.2.2, relative to freestream 
(inertial) reference datum. Translational displacements (X, Y, 2)^ of the 

origin of the fuselage coordinate, and rotational displacements { $, 6, •*}„ 


about the fuselage axes describe these degrees of freedom. See Figure 8. 

As mentioned in Section *.2.7, the freestream set may instantaneously assume 

any reference point; therefore, only the time derivatives of (X, Y, 2) and 

Oh 

( *, 6, ♦) have significance. In order to locate the direction of the 
n 


gravity vector relative to the hub, a running calculation of the Euler 

angles 6-, i>- must be made. Since these are not degrees of freedom and 
t r* i. 


therefore not calculated in the equations of motion, they must be calculated 
outside the dynamic equations as the tine history proceeds. When the initial 
orientation of the hub is defined, 9^, and 4>£ are known and their 


changing values nay be calculated by integrating the hub rotation rates in 
the earth or freestream axes. 


h.2 .2 Rotor . - The rotation for the rotor degree of freedom is defined 

as motion of the rotor coordinate system relative to the hub axis system. 
This is shown in Figure 6. This figure also indicates the change from 2 
down to 2 up axis, which is equivalent to a 180 -degree positive rotation 
about the Y axis. Note: Rotor rotation also includes blade feathering 

from swashplate rotation in addition to blade root rotation. 


k.3.3 Shaft or transmission deflections . - Shaft or transmission degrees- 
of-freedom are defined as motions of the hub coordinates relative to the 
shaft axis system. Hence, as shown in Figure 4, hub notions are dependent 
variables which are functions of the shaft deflections. 



ITEM 


SYMBOL 


TYPE OF MOTION 


FUSELAGE 
PRINCIPAL AXIS 


Xgp, y of- Zof TRANSLATION AND ROTATION 

& a WITH RESPECT TO 

F*YF*'lF INERTIAL REFERENCE 


ROTOR 


ROTATION OF ROTOR SET 
WITH RESPECT TO HUB 
AXES 


SHAFT OR TRANSMISSION 
DEFLECTION 


X S Y* 2 * 
< V°S-'-'S 


DEFLECTION OF HUB SET 
WITH RESPECT TO SHAFT AXES 


BLADES 

(n = 7 MAXIMUM) 


A 1«t* A 2n- A 3n 
PPHn 


BLADE RENDING MODES AND 
FEATHER/PITCH HORN 
BENDING OR TORSION WITH 
RESPECT TO BLADE 
ROOT AXES 


SWASHPLATE 


< *’SP- e SP 

Z SP 


SWASHPLATE AXES MOTIONS 
WITH RESPECT TO HUB 
AXES 


Figure li. 


Tegreen of 


freedom. 
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^.3.** Blades. - Each blade's motion relative to the rotor coordinate system 
is defined in terms of four generalized coordinates. These consist of three 
blade bending modes and a combined feathering, pitch arm bending mode, or a 
torsion mode. 

U. 3.1.1 Blade bending . - Blade motion due to blade bending is defined by 

the generalized medal coordinates A which typically represent a coupled 

mn 

first deplane bending node, a coupled first flapvise bending mode, and a 
coupled second flapvise bending mode. Ordinarily in a modal analysis, the 
effects of centrifugal and structural stiffness are lumped together into a 
generalized stiffness which is simply the modal natural frequency squared 
times the generalized mass. In contrast to this, the REXOR II analysis sep- 
arately treats the strain energy or structural stiffness in each mode and 
the stiffening due to the centrifugal force field. This provides the capa- 
bility of being able to account for the periodic variation of stiffness in 
the modes due to the reorientation of the centrifugal force field with re- 
spect to the blade principal axis due to variations in blade angle. This 
feature can be important in the study of subharmonic stability where the 
periodic variation of coefficients may be important, but it also permits 
being able to make rather large changes in rotor speed and collective blade 
angle without having to change blade modal data. 

Mode shapes and natural frequencies are initially determined for a twisted 
blade at or near the collective blade angle and rotor speed to be analyzed. 
Such effects as precone, blade sweep, blade droop, and blade angle varia- 
tion are included in the REXCR II analysis and couple the initially' orthogonal 
modes. The elastic bending contribution due to the modal deflections is 
calculated relative to the blade's static shape. 

As previously noted, the blade modes are initially defined at some refer- 
ence feathering angle, ♦ As time progresses in the analysis, the blade 

feather angle varies about this reference position. The mode shapes are 
correspondingly transformed to account for the difference between the in- 
stantaneous feathering angle and the reference feathering angle, at the 
same time accounting for other effects such as the static and instantaneous 
shape of the blades. This yields the modal coefficients (partial deriva- 
tives) that relate blade element motion to the blade bending generalized 
coordinates as a function of time. 

The vertical and inplane blade element variational motions, 6Y. and 6Z., 

1 l 

can be written as follows: 


3Y. 


3Y. 


3Y. 

l 


<Y in = 3*r ( V t)5A ln * 3aT ; V t,<4 2 „ * 5a 7 ( V t)SA 

In ?.n 3n 


13) 


3n 
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and 


6 Z. = t 

in 3. 


az )z 

( v t)6A m + * 


-IP 


2n 


3Z 

■' t)SA 2., *• IaT ( V t,SA 3a 

3n 


(*») 


where the given cr input partial derivatives are the true nodal coefficients 
of the orthogonal modes for the blade in an undeformed shape, with no static 
geometry accounted for, and at the rotor speed and collective angle for which 
the blade modes were initially calculated. 

The orthogonal bending inodes used in the analysis are illustrated in Fig- 
ures 12, 13, and lU. Observe that the root boundary conditions for the 
modes may be cantilevered or articulated. 

Hote that in addition to the normal bending responses, Y^ and Z^, the span- 

vise motion of each blade element is also determined, and blade feathering 
due to pitch-lag and pitch-flap kinematic coupl ng effects are also accounted 
for in each blade bending mode. This feathering is added to that due to 
svashplate motion as is blade feathering due to flexibility. 

This modal data is developed to the form used in the blade equations in 
Section U.5.5. The discussion of modes is carried on from a math view- 
point in Section 5- 1 *- 

U.3.U.2 Pitch horn bending - dynamic torsion . - The remaining mode per 

blade, pitch horn bending, is comprised of either a blade feathering drive 

flexibility with a torsionally rigid blade cr an uncoupled torsion mode. 

Examining the first alternative, the svashplate position determines the 

primary blade feathering motion. In addition, the linkage between the 

svashplate and the blade (see Figure 15) has flexibility in the pitch link, 

pitch horn, and cuff. The feathering or pitch horn bending degree-of- 

freedom therefore can be rigid blade feathering motion outboard of the 

blade cuff coupled with a net inboard stiffness. Inboard of the blade 

cuff, feathering flexibility results from the pitch link, pitch link 

bearings, pitch hern, and cuff. The relationship between blade feathering, 

, and motion of this degree-of-freedom, , is defined as the partial 
rn rrH 

derivative. 


M 

36 PH /n 


37 










can be a distributed torsional 


Alternatively, this degree of freedom, 8 pii 

response of the blade based upon defining an uncoupled dynamic torsion mode. 
The selection of the degree-of-freedom representation is made on the basis 
of the type of analysis being performed. The mode defined is uncoupled in 
the sense that it is not a function of the flapping or lead-lag modes. 

An optional quasi-steady torsional response of the blade may be used in 
conjunction vith pitch horn bending. This is superimposed on the rigid 
blade feathering and permits a distributed torsional response alternative 
of the blade reacting the spanwise variation of applied torsional momen s 
from aerodynamics, coriolis, and centrifugal force terms. The blade tor- 
sional response at the ith blade station is computed from the following 
equat ion : 


Ti 


i— P 1 

t T S+1 J root 




M (x)dx 
$ 


(5) 


where S is the Laplace operator, and is the time constant associated 

± 

with blade torsional response. This equation is implemented numerically 
in the REXOR II program. 

To aid in program trouble shooting the pitch horn bending representation 
(with or without quasi-static torsion) may also be operated as a qui si- 
static degree of freedom without second -order response. 

-.3.5 Swashplate. - The swashplate has three degrees of freedom: 

' bi 

6 , and Z . Rotations $ and 6 are Euler angles defining the orienta- 

bi bi OJr bi 

tion of swashplate coordinates relative to the hub. Likewise, the transla- 
tion Zgp defines vertical displacement of the swashplate relative to the hub 

axis. These nre shown schematically in Figure 10. 


U.k General Motion and Coordinate Transformations 

In development of the equations of motion, it is convenient to write the 
forces, moments, velocities, and accelerations in coordinate systems 
related to separate elements of the system. Consider the concept of 
general space motion of a particle. 



U.U.l General case of space notion . - For the general case of space motion, 
a particle, p, moves vith respect to a reference axis system which is, in 
turn, in motion vith respect to a fixed coordinate system. This is illus- 
trated in Figure 16 where the fixed or inertial coordinate system is des- 
ignated by capital letters X, Y, Z, and the moving coordinate system is 
designated by lower case letters x, y, z. The moving coordinate system is 
rotating at an angular velocity, w. The vector u may, in general, vary in 
magnitude and direction, both of vhich can be referenced with respect to 
the fixed X, Y, Z axes. 

Thus, the absolute motion of the particle p, referred to the inertial 
X, Y, Z axes, is equal to the motion of the particle relative to the 
moving coordinate axes x, y, z plus the motion of the moving axis system 
with respect to inertial space. 


u) . 



Figure l6. - General case of space motion in terms of moving coordinate 
axes x, y, z and inertial axes X, Y, Z, 


Ul 



To vis ual ize the motion of the particle p, let its motion with respect 
tc the '■'.oving axis system be indicated along a curve s fixed in the 
moving axis system, x, y, z. An observer sitting on the moving axis 
system would therefore see only the motion of p along the curve r.. 

From Figure l6, the position of p relative to the x, y, z axes is represented 
by the vector 


r R = xi + yj + zk 


( 6 ) 


where i, J, and k are unit vectors along x, y, z, and therefore must 
be treated as variables due to their changing direction. Differentiating 
-► 

r R results in 


— f f T t -r *r 

• • • di d.1 dk 

r = xi + yj + zk + x-n- + y^ + z — 


dt 


dt dt 


(?) 


di -► -► d 1 dk 

Since — = w x i, -rf = w x j and — = m x k, this expression can be 
at at at 


written as 


• • • • . . 
r = xi + yj + zk + u> x (xi + yj + zk) 


( 8 ) 


or 


r = r + u x l 


(9) 


In this equation, the first term, r represents the velocity p relative 
to the rotating axis, x, y, z. The second term, w x r, is the velocity 
of the point in the moving coordinate system due to the rotation w. The 

absolute or inertial velocity R of the point p is obtained by adding 

■+ 

the velocity of the origin of the moving axis system to or: 


R = RQ + r + u)^.r 


where 


w = pi + qj + rk 


and 

1*2 


uj = pi + qj + rk 


( 10 ) 



The inertial accelerations of the point p can now be determined by 
simply differentiating this expression with respect to time. Performing 
this differentiation yields 


R 


R 0 + r + a)xuxr + a)xr + 2u)xr 


( 11 ) 


where the terms a> x w x r and Jxr represent accelerations of the 

-*• 

coincident point in the moving axis system, r is the acceleration of p 
relative to the moving axes, x, y, z, and 2<u x r is the coriolis 
acceleration which is directed normal to the plane containing the vectors 

u and the relative velocity r, as given by the right-hand rule. 

The vectors expressed in the preceding equations are in the most general 
form for defining the motion of a particle moving in a moving coordinate 
system. All special cases can be deduced T'orn these equations. 


For convenience, the time derivative equations can be expanded in matrix 
form. The inertial or absolute velocity and accelerations of the particle 
p, written in expanded matrix form, are given by: 
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Performing the indicated matrix multiplication gives: 
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U.U.2 Coordinate transformations - Euler angles . - To describe motions in 
one coordinate system in terms of motions in another coordinate system, 

Euler angles $, 0, and 4 > with the appropriate subscripts are introduced. 
These angles can be applied to define the rotation of one coordinate system, 
x, y, z, relative to another coordinate reference frame, X, Y, Z. Since the 
development contained in this report utilizes these angles in relating 
coordinate systems, a brief explanation is given here. 

Rotational displacement of a coordinate system can be represented by the 
three rotational displacements <p , 0, and ip, as shown in Figure IT . The 
order of rotation is not important as long as the sequence selected re- 
mains consistent and the reverse order is used when rotating back to the 




original position. In this analysis, the rotations start with displace- 
ment ♦ about the x axis, then a rotation 8 about the nev y axis, 
followed by a rotation $ about the new or final z axis unless geometry 
or physical considerations of the modeled part dictates another order. 

This Mans the (X, Y, Z) coordinates can be rotated into the 

(X, Y, Z). axis system as follows: 
b 



and the inverse transformation can be written as 
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Sty inspection, then, it can be seen that 

inverse of [f] = transpose of [tJ 


( 20 ) 


or 


[vj 1 ■ M - [vj 


( 21 ) 


Carrying out the indicated matrix multiplication yields the transformation 
matrix [T] : 
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Using this transformation, the inertial velocities and accelerations of a 
point or particle be written in one coordinate system in terms of those in 
the other coordinate system as follows: 
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and: 


and inversely. 


and 
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i,.U„3 Angular velocities and accelerations - general . - ror the general 
case, consider the coordinates in the previous section, and let (p, q, r) 

a 

and (p, q, r)^ be the respective angular velocities of and about the 

(x, y, 2 ) and (s, y, a) axis systems. Also, assume that the Euler angles 
& 0 

are varying with time ($, 0, and ♦), and let (x, y, 2 ) be the reference 
coordinate set with (x, y, 2 )^ coordinate set moving relative to it. This 
is illustrated in Figure 18. 
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From this figure, the following can be written. 



Differentiating tnis expression with respect to time results in angular 

accelerations (p, q, r). in terms of the reference coordinate system 

b 

angular velocities and accelerations. This results in the following: 
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These equations represent a general form for defining angular velocities 
and accelerations of one axis system rotating relative to another axis 
system, which in turn is in motion. 

A special case is the angular velocities of system b with zero Euler angles. 
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k.5 Relative Motions and Transformations Used in the Equations of Motion 

In this section the inertial linear and angular velocities and accelera- 
tions of major components of the vehicle, are presented. Also included is 
the development of coordinate transformations that relate motion in one 
axis system to another. Motion of the principal reference axis system in 
relation to the earth is described. Motion of each component or reference 
axis system is then defined in terms of the degrees of freedom. 

U.5.1 Fuselage motion in inert ial space . - At each instant in time the 
fuselage axis (Section 2+. 2.1) is related to an inert ial coordinate axis 
system. Inertial accelerations of the fuselage axis system are defined 
by the vector 


r 


0 

v. J 


(29) 


where the quantities represent the total inertial acceleration of the 
generalized coordinates of the vehicle as defined by motion of the 
principal coordinate axis system. 
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Orientation of this system relative to the earth is specified hy Euler 
angles # , 0 , and as seen in Figure 8. The sequence and def ini- 
tion of these angles is ♦ (yaw), 0 £ (pitch), # £ (roll). Hote that the 
sequence of rotations is opposite to that given by Figure 17 • The 
angular rates, p p , q p , r p , of the fuselage or principal reference axis 
system with respect to the inertial coordinate system can be written as 
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This equation can be rewritten to solve for # , 0 , and as 
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The Euler a;..:es defining orientation of the principal reference axis 
system with respect to the earth is next obtained by integrating the rates 
with respect to time, or 
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Angular velocities of the fuselage, with respect to the inertial axes 


reference system, p p , q p , r p , are defined in terms of the degrees of freedom 


as 



dt 
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P 
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r t « 0,F 


( 35 ) 


where 



( 36 ) 


Linear velocities of the fuselage or principal axis system are now determined. 
The first three quantities of the fuselage axis acceleration vector represent 
the linear interial acceleration of the fuselage. For a system in motion, the 

inertial acceleration, a^, at the origin of the system is defined, 

based on the vector algebra of Section U.U.l, as 



a„ + u x V, 


( 37 ) 


-*■ dV_ 

** 0 

where a Q is is the rate of change of velocity, V Q , of the origin 

of the moving coordinate system and w is the rotational velocity of the 
moving coordinate system, both relative to the earth. Now defining: 
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gives 
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From this equation, then, the rate of change of velocity of the moving 
coordinate reference system becomes 
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This set of accelerations and the time integral represent airflow accelera- 
tion and velocity incident on the helicopter. 


a separate set of accelerations is carried through the analysis which contain 
the acceleration due to gravity. Ordinarily, gravity is treated as a force 
of mg on the right-hand side of the equations. However, the gravitational 
term can be accounted for by defining 



(M) 


where g^, g^, g^ are the three components of the gravity vector to be 
defined. The acceleration on the left may be defined as being in earth- 
inertial, El, axes. 
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The logic behind this substitution is as follows. For a rigid body in 
motion, the equilibrium equations can be written as 


mX = m(u + qw - rv) * F 
mY = m(v + ru - pw) « 
mZ = m(w + pv - qu) = F„ 

LA 

where 
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( 1 * 3 ) 


F , F v , and F„ represent the external forces acting on the body, exclusive 
X i £ 

of gravitational forces. 


Subtracting the gravitational vector from each side of the previous 
equations yields: 


m(X - 
m(Y - gy) 
m(Z - g z ) 
which by inspection gives 

X - 
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m(u + qw - rv - g^) = ? x 
m( v + ru - pw - g y ) = F y 
m(w + pv - qu - g^) = F z 

= u + qw - rv - g x 
g y = v + ru - pw - g y 
g 2 =w+pv-qu-q z 
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( 1 * 5 ) 

(1*6) 
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( 1 * 9 ) 
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Rearranging these equations yields: 


u = (X - g x ) - qv + rv + g x 


- (Y - g v ) - ru + pv + g 


w = (Z - g z ) - pv + qu + g 2 


(50) 

(51) 

(52) 


The first terms on the right side of the equation are identified with the 
proposed gravitational acceleration definition of Equation 41. 

Making the substitution: 
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In this equation, the accelerations X^, Y^J, 


-Ei 

and Z_„ are the degree- 

Or 


of-freedom accelerations of the principal reference axis system used in 
the REXOR II analysis. These accelerations represent the inertial accel- 
erations plus the equivalent accelerations of the reaction force to 
gravity. Thus, gravity is an equivalent acceleration applied to the 
reference coordinate axis system. Via coordinate system referencing, 
every mass element on the vehicle is therefore acted upon by this accel- 
eration. This avoids including gravitational force as an external force 
individually applied to each mass element. 

The gravitational vector at the fuselage is simply the gravity vector in earth 
axis transformed to the fuselage axis system through the Euler angle rotations 

<f> E » 0 £ . and <J> E . Or 
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where 



0 

cos$ E 

-sin$ E 


0 

sin *E 

COS(> E 


COS0 £ 

0 

sin0 E 


0 

1 

0 


-sine E 

0 

cos 6_ 
E 


COS^/g 

-sin* E 

0 


sinijig 0 

cosing 0 

0 1 


( 55 ) 


The velocities of the principal axis system are obtained by integrating 
the rates of change of velocity with time, or 
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These velocities in earth coordinates can be written as 
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which can be integrated to give the position of the system relative to the 
earth. Doing this yields 
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U.5*2 Hub motions in inertial space . - The geometry of the hub, shaft and 
fuselage is described in Sections U.2.1 through h.2.3. Forming the motions 
of the hub first requires knowledge of the shaft set motions and shaft 
generalized coordinates. The transform from fuselage to shaft set,|T F , 

involves rotor tilt alignment data $ and 8 . Elastic motions of the 

o o 

s s 

transmission suspension, and consequently the hub, are described by the 
generalized coordinates (X, Y, 2, $, 0 , s which are measured with respect 

to the shaft (S) set. The transform from shaft to hub, J? s jjj is a function 

of the generalized coordinate angles ( <f>, 6, <iOg. 

The development starts with fuselage to shaft set relations: 



Noting the following transformations: 
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Applied to hut set equations: 
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1.5.3 Motion of rotor coordinate axis . - The rotor coordinate axis system 
is shown in Figure 6. Note that the rotor coordinate axis system is rotated 
180 degrees about the Y axis relative to the hub axis system at the time 
when the rotor is at azimuth position zero. That is, X and Z change direc- 
tions. The rotor coordinate system then rotates through the angle \j> R from 
this position. 

The sequence of rotation in going from hub to rotor coordinates consists of 
first a 180-degree 0 rotation, followed by the * R rotation. Following the 

convention established in Section t.U.2 for Euler angles: 
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where the 

last matrix represents the 180 -degree 

0 

rotation. 

The next 




ft 







matrix is 

the rotor rotation. 


s 

dt. 






Since the origins of the rotor coordinate system and the principal refer- 
ence axis system are coincident, the linear velocities and accelerations 
of the origin of the rotor coordinate system can be directly written as: 
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and 



( 73 ) 


Noting gravity has been treated as an equivalent acceleration in the hub 
generalized coordinate accelerations. This sane equivalent acceleration 
El 

is included in (X„, Y_, Z-) D , the rotor coordinate accelerations. 

Q 0 u n 


The angular velocities, p R , 

the rotor coordinate system are determined; again noting the rotation 
order. The rotor coordinate system angular velocities are: 


q R , r R , and accelerations, P R » 4 R> r R of 
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likewise, accelerations of the rotor coordinate system are: 
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The above equations then define the coordinate transformation from hub to 
rotor coordinates; and rotor axis system linear and angular velocities and 
accelerations in terms of velocities and accelerations of hub and the rotor 

degrees of freedom t . 

K 
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U. 5 . 1 * Blade coordinate relative to rotor coordinates . - Since each blade 
has its own blade reference system, as shown in Figure 6, the and Yg^ 

and Yg axes azimuthal ly by an angle 

defined by the equation 

_ -2w(n-l) (76) 

? BLn b 

where b is the number of blades and n is the blade number. This equation 

and the Y R axes are coincident. 

The transformations between the rotor coordinate axis system and the blade 
coordinate axis systems are defined by the equation 


states that the Xg^ and X^, and the 


axes are rotated with respect to the Xg 
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Note that these equations define blaae one as being straight aft at time 
zero. 
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In the blade reference axes, the velocities and accelerations of the 
origin of the blade reference axis system become: 
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Likewise, the angular velocities and accelerations of the blade reference 
axis systems become: 
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and 
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it. 5-5 Blade element motion . - The following blade motion description, due 
to the involved nature c' the geometry, is rather lengthy. First, in this 
development, the motion of the blade with respect to the relative blade 
coordinates is given. This motion is the sum of static and modal deflec- 
tions. Then the relation to free at ream coordinates is computed. Partial 
derivatives sire extracted from the transformations for use in the equations 
of motion of the blade in Section 5-6. 

The blade element motions for the n th blade are defined relative to the 
blade (BLn) coordinate reference axes (Figure 6). The blade element rela- 
tive motions are functions of the static shape, of blade feathering and 
torsional deflection, and of blade bending of the coupled inplane and 
flapping modes. 

The static shape includes such items as blade twist, hub precone 

angle, 8„, blade droop angle relative to the precone angle, v, blade 
sweep angle, t^, feathering axis precone, 8p^, the blade feathering 
angle, and the blade elaaent center of gravity location. 

The blade motions about this static shape include the effects of the 
three blade bending modes, A^, and A^, blade feathering, 
and blade torsional deflection, 4 ^. 

The blade elaaent motions are now defined. The blade static position in 
the blade reference axis system is first developed. The blade bending and 
feathering deflections are then introduced. Both deflections and slopes 
are developed and then these equations are differentiated with respect to 
time to obtain the bxade elaaent linear and angular velocities and 
accelerations. 


The blade element linear motions are developed in blade (BLn) coordinates 
and the blade elaaent angular velocities and accelerations are developed 
in blade elaaent (BLE) coordinates. The coordinate transformation 


matrix 


[ T BLn-BUs] 


‘BLn BLEJ * s a ^ so to peraj t the transformation of the 

inertial velocities and accelerations from one axis system to the other. 
The development of the blade relative motion equations now starts with 
the description of the shape of the blade. 
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U.5.5-1 Blade static shape. - Blade elemental motion is defined as motion of 
the blade element reference axis system which has its origin at the blade 
element center of gravity. The blade aerodynamic reference axis is selected 
as the 1/U chord. Likewise, the geometry and dynamics are referenced to the 
1/U chord, though any reference line could have been used. Starting with the 
straight untwisted blade with the blade 1/1* chord lying along the axis 

as in Figure 19, the blade element eg and blade element coordinate axis system 
origin are defined by the coincident point defined by the vector 
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in blade coordinates. 


Thd dimension X__(i) nT is the undeformed spanwise 

Lli bLn 


location of the cg/blade element origin. The dimension Y„ 0 (i)_ T is the 

lu tilM 

chordwise location of the c.g. /blade element axis system origin forward of 
the blade lA chord and Z^i ) BLn is any vertical offset of the c.g. /blade 

element origin with respect to the reference chord plane of the blade. 


Now, introducing blade twist by rotating about the X RTn axis through the 
local blade twist angles. Figure 20, results in: 
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The Roman numeral subscript I denotes the first of a sequence of static line 
transformations . 
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At this point the subscripting, BLE will be dropped to simplify the develop- 
ment. Rewriting the above equation, we have: 
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Introducing blade coning, 6^, results in the location of the blade as shown 
in Figure 21. This results in: 
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The next item of static geometry that is considered is blade droop, y, and 
then blade sweep, t . These rotations are shown in Figure 22. Note that 

since the blade sweep and droop angles are introduced at a distance out 

on the blade, it is first necessary to transfer axes to this location before 
making the rotations. Therefore, the blade displacements outboard of 
Station X_„ become: 
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At this same station, provisions are introduced to allow for offsets of the 
blade in both the vertical and horizontal directions by Zj 0 g ^jog* 

respectively. These offsets are shown in Figure 23* These offsets represent 
displacement of the blade 1/U chord with respect to the blade precone line at 
blade station Xgy 


Introduction these offsets, then, and transferring back to the center of rota- 
tion through X„„ results in the description of the blade displacements out- 
i>W 

board of station including the effects of precone, sweep, droop, and 

offset of the blade from the precone line. 
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and Z 


jog 


Figure 23. - Introduction of blade 1/r chord offset, Y. 

with respect to precone line. J K 
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At this point, a reminder that the prior development represents the blade 
displacement inboard of Station X~ w and the above equation outboard of 

Station Xg y . Therefore, inboard of Station Xg^: 
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With this in mind, the remaining developing of including the effects of 
feathering axis static precone and blade reference feather angle in describing 
the static blade position continues. No distinction will be made in the fol- 
lowing developments between inboard of Station X gw and outboard of Station X^ w 

Figure 2l shows how blade feathering is introduced. The axis system is trans- 
lated to a point p which is located at the intersection of the precone line 
and the feathering axis. The location of this point is a distance £ along 

the cone line, as shown in this figure. The blade is first rotated to the 
feather axis; then rotated about the reference feathering angle, <t> , the 

feathering angle for which the blade modes are defined. Doing this results in 
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FEATHERING AXIS 



Figure 2b. - Point p and feathering axis precor? (3^. 


This equation defines the looatior of the static shape of the blade in an 
axis system with the y-axis horizontal and the x-axis aligned with the blade 
static feathering axis. Transforming now back through the feathering axis 
precone angle and translating back to the rotor shaft centerline results in 
the static shape of the blade defined in blade coordinates, or 
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where subscript S refers to blade static or undefcrmed shape. Combining 
equations developed so far results, then, in the following two equations which 
represent the static shape of the blade for both inboard and outboard of 
blade station X„ TT . 
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These two equations then define completely the static shape of the ' lade. 

The dco'lopment will not proceed to include the blade bending or elastic de- 
formation. .iowever, before proceeding with this, the static location of the 
blade feathering bearing is defined since these will be used in the develop- 
ment that follows. 
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Referring to Figure 25, it can be seen that the static position of the inboard 
feather bearing location can be written as: 
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The static location of the outboard feather bearing is: 
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With these definitions, the analysis will rot proceed to include the effects 
of blade bending, blade feathering, and torsional deflection. 


ZsLe 



Figure 25. - Static feather bearing geometry. 
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4.5.5. 2 Blade shape - elastic deformation . - In the foregoing development, 
the analysis has proceeded in a completely rigorous fashion. At this point, 
though, a departure from a completely rigorous simulation of the elemental 
blade motions will be made. It will be assumed, as far as blade elastic de- 
formation is concerned, that the cosine of angles, like precone less droop, 
blade sweep, elastic flapping, and elastic inplane slopes, but not blade 
feathering is approximately equal to 1, and therefore, the blade elastic de- 
flections, y and z, in blade coordinates, will be assumed to be equal to those 
in the static blade element coordinates. This assumption in a reasonably 
valid assumption and is completely consistent with standard practice in the 
mathematical representation of blade element motions. 

Additionally, as far as the effect on structural axis reorientation due to 
blade 9 rotation, the effect due to blade leastic twist is considered to be 
small compared to that due to blade cyclic and collective feathering. Also it 
will be assumed that the contributions to blade Y and Z motion are small due 
to blade torsional motion, other than that due to local center of gravity 
offset. 

With these as." motions in mind, blade elastic trading will now be introduced. 
The contribution to elastic blade bending is si-^ly 


0 


1 

o 

o 

o 

1 


r 1 
A ln 

f BEUD 

► — 

Y Y Y 

1 2 3 

< 

A 2n 

"BEND 

w A 

BLn 

7 *7 7 

"1 ^2 Z 3 


A, 

3n 

» 4 


( 96 ) 


Note that X or spanwise motions are not included in this equation. Blade 
spanwise motion will be determined separately by utilizing blade slope data 
to determine the change ir. the projected blade length upon the blade X axis. 
With this in mind, the total Y and Z blade motions including blade bending, 
but not yet including blade feathering or blade elastic twist , is strictly 
the sum of the previous static line expressions and the modal deflection. 
Blade torsional aeflectior is treated as an independent degree of freedc-n, 
and therefore is not included as part of these blade modes. Combining the 
previous static deflection with the modal deflections gives: 
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^.5.5.3 Blade feathering . - Blade feathering is relative to the reference 
feathering angle The feather angle, then, as far as blade motion is 

Khr 

concerned, is due to the difference in the total feather angle and the 
reference feether angle 

The blade feathering motion is introduced similarly to the way the blade 
reference feathering angle was introduced, except that the feather axis slopes 
are due to the static position as well as due to elastic deformation in both 
the flapwise and inplane deflection. 

If we let and Y 'fa re P rese,rt the instantaneous vertical and inplane slopes 

of the feathering axis, then transferring to the inboard feathering bearing, 
making the rotations through Z'^ and *'fa to the feat ^ erin € axis, rotating 

through the delta feather angle - ($p -t p£p) of-A^, rotating back through 

-Y* and -Z'_ ,, and then transferring back to the BLn axis system results 
rA rA 

in the definition of the displacements in blade axis coordinates. 

However, before proceeding with this, the feathering axis slopes and 

are defined. The slopes are simply defined as the difference in the 

total static and elastic deflection of the outboard and inboard feather bear- 
ings divided by the spanwise distance between the bearings. Then from Fig- 
ure 25 and the bearing static location equation: 
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where in terms of the static anl modal deflections 


and 
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In the development that followr, the time derivatives of Y* ^ and Z'p A are 

required, so therefore, they are now defined. Taking the first and second 
time derivatives of the slope equations yields 
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where 





(109) 


Transferring the blade displacements as indicated above to the inboard feather 
bearing, transforming to the feathering axis, and performing the feathering 
rotation as discussed earlier, yields the following equation which defines 
the blade displacements in blade axis coordinates: 
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This equation then gives the blade displacement in blade coordinates, includ- 
ing the effects of the static shape, blade bending, and blade static twist. 

The effect of blade elastic twist is now considered. 

4. 5. 5. 4 Blade elastic twist . - Blade motion due to blade elastic twist is 
accounted for by going back to the static twist equation. Blade elastic twist, 
$ is assumed to be directly superpositionable with blade static or blade 

pretwist, except that the static pretwist takes place about the 1/4 chord, 

and the blade elastic twist takes place about the blade element shear center. 
This is shown in Figure 26. From this figure it can be seen that previous 
static twist equation car. be rewritten as: 
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a) BLADE PRETWIST. <>M m . ABOUT BLADE REFERENCE AXIS 



b) BLADE ELASTIC TWIST. *<0 T ABOUT BLADE SHEAR CENTER 

Figure 26 . - Blade static pretvii, and elastic twist, 
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I+.5.5.5 Final blade element Y, Z displacement equation. - Suostituting the 
above equation in the previous development sequence yields the blade dis- 
placement equation which includes the effect of the static shape of blade 
bending, of blade feathering, and of blade elastic twist. 

However, beffe proceeding with these substitutions, the following column 
vector is defined to simplify the notation. 
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The total blade element displacement equation becomes: 



where: 
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Note that for convenience of using the condensed matrix notation discussed 
above, the most general vectors for such terms as £^, Xg^, and have 

been used. As can be seen in this equation, these have all been treated as 
full vectors. Making the appropriate substitutions of course will result in 
the expressions previously obtained. 

It is noted that the equation is written for the relative displacement of 
points on the blade outboard of Station Xg W - Inboard of that station, the 

displacements are determined from the previous onboard equation or simply by 
zeroing out such terms as | r iog| j r < 3 vfj substituting unit diagonal 

transformations for JT t J and ft] in the full equation. Following either 

approach yields the blade displacement equation for points inboard of 
Station X-.,; or 
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The ith station blade displacements, Y and Z, in blade coordinates for points 
on the blade both outboard ara inboard of station X ^ are then defined. 

1.5*5 .6 Blade element Y ana Z relative velocities and accelerations . - The 
blade element coordinate axis system linear Y and Z velocities relative to the 
blade reference axis system can be found by differentiating the position 
equation with respect to time. Note that no distinction will be made at this 
point between outboard or inboard of station Xg W , but using the equation for 

displacements outboard of this station and as discussed earlier, zeroing out 
certain terms, results in the equations for velocities or accelerations of 
points inboard of that station. 
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Note in the above equation that the matrices are not time derivatives of 

the £T;.J matrices but are derivatives of the transformation matrices with re- 
spect to the transformation angle £. . This is arrived at by making the sub- 
stitution that: 


and 
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Taking the time derivative again of Equation 113 yields the blade element 
Y and Z linear accelerations relative to the blade reference axis systc-a. 
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These equations define the blade element relative displacement velocities and 
accelerations, respectively, required by the blade inertial velocity equations 
developed shortly. Rote that in the preceding development these ^equations are 

r 3r 

written for the nth blade, and with the exception of the 

[%]• [ T >} N- NJ- M’ l rsw l' l rjo8 l' N 

ces, the terms are all blade dependent. Remember, also, that inboard of 



matrices, the terms are i 

M tte [\] " d M 


matrices are unit diagonal. 


U.5.5.7 Blade element slopes . - The blade element Y* and Z* slopes are 
determined by differentiating the deflection equation with respect to the nth 
blade radial distance, 1^. These formulations are used for quasi-static 

torsion formulation and output. Performing the required differentiation for 
points along the blade reference line: 
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2*. 5. 5. 8 Transformation from blade-to-blade element coordinates . - In this 
section, the transformation matrix for the blade root to the ith blade element 
motion will be developed. Each blade will have its own transformation matrix 
for each ith blade station. The transformation matrix will initially. be 
developed as the transform from blade element to blade coordinates, FuLE-BLn I ’ 


The transformation matrix. 


[ T BLE-BLn] 


, can be de 'eloped by referring to the 


development of the deflection equations. The fi’-st rotation from blade ele- 
ment to blade coordinates is through the combined twist angle, - 4 ^; the second 

rotation is through the negative of the precone, 8^; the third through the 
negative of thw sweep and droop angles, Tq, y ; the fourth through the feather- 
ing axis angle, 6p^; the fifth through the negative of the reference feather^ 
ing angles and the sixth back through the negative of the feathering 

axis p-econe angle, 8 FA « 

These rotations then define the transformation from blade element to blade 
coordinates, including the effects of the static shape of the blade, pretwist, 
precone, sweep, droop, etc. Also included is the effect of blade elastic 
twist. Again note that for stations inboard of Station the sweep and 

droop angles, and y. respectively must be set to zero in the formulation of 

the transformation matrix as in the definition of the blade displacements and 
blade slopes. This portion of the transformation matrix which includes the 
static blade shape and combined twist is defined as follows: 

[w.].-[M [vr [*,.]] [[\]’ t-r M T HI 

(12U) 


The next two rotations from blade element to blade coordinates are due to the 


elastic blade bending slopes. 


Since y ’bend and 


BEND 


are motions of the 


blade elements with respect to the blade, then to transform from blade element 

to blade coordinates requires negative rotations of Y' and Z* __ to be 

cLIvD BEND 

included. Finally, the blade feathering rotation from the reference feather 
angle must be included. The final transformation then, from blade element to 
blade coordinates, is defined by the following equation: 
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Also: 
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The inverse or transpose of this equation yields the transformation from blade 
to blade element coordinates, or: 



MMHMW WM 



bb-d] [ Tz Va] T W 
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inboard of station X„ tr . 
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U. 5 . 5.9 Blade element angular velocities and accelerations . - From the fore- 
going discussion, the blade element angular velocity vector can be determined. 
Starting with the angular velocities (p, q, r ) BLn of the blade reference axis 

system and systematically and progressively transforming these velocities 
through each axis rotation and adding the respective angular velocity asso- 
ciated with each of the indicated angular rotations, results in the following 
equation for the blade element angular velocities. 
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Bote that in this equation, starting on the right-hand side vith the quantities 
in the innermost brackets, the blade reference system angular velocities are 
first transformed through the increment of feathering axis flapping slope due 
to bending, 2 ' and then the feathering axis flapping angular velocity, 

is added. Minus is used since 2' is a negative 0 rotation. Next, the 

resultant W vector is transformed through Y' and Y' is added. This is 

fA rA 

then transformed through the delta feathering angle, A$p, ^ the feathering 
angular velocity, ip, is added. This is then transformed back through the 
increments of feathering axis slopes due to blade bending, giving the vector: 


9U 




which represents the blade element angular velocities due to combined blade 
feathering and blade reference axis system angular velocities. 


Next, the effects of blade bending r.t each blade station are introduced. The 
above vector is first transformed through the local blade element flapwise 
bending slope, Z ' BEND » and then the angular velocity, ” z 'bend* * s a< * de< *‘ This 

result is transformed through the blade element inplane bending slope, Y* ___, 

BEND 

and the inplane angular velocity due to blade bending, is added, re- 

sulting in the total vector less the initial transformation string. This 
vector represents the blade element angular velocities due to the combined 
effects of the blade reference axis system angular velocities of the blade 
feathering angle and of the blade angular velocities due to blade elastic 
bending. The remaining transformations then include the static effects of 
the blade feathering axis precone, Bp^, the blade reference feathering angle, 

it'ppp, blade sweep, blade droop, y, and blade or hub precone, 8 qi and the 
combined effect of blade static and elastic twist, represented by Finally, 

the blade elastic twist angular velocity, <j>_, is added, giving the total 

J1 

blade element angular velocities, J q v 


Also note, as indicated before, the matrix jT^J T t ' has the value calculated 

if X is greater than and has the value of unity if X is inboard of 

station X . 
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where: 
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and: 
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The blade element angular accelerations can now be determined by differentiat- 
ing this equation with respect to time. Again, as in the case of the angular 
velocities, the contributions due to time derivatives of the feathering axis 
flapping and inplane slope changes due to bending are neglected. With this 
assumption, the time derivative is: 
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k. 5 .5.10 Blade element X potions . - In the previous development, the equations 
did not account for the blade element displacement, velocity, and acceleration 
in the spanvise or X-direction. The method used to define these is one of 
taking the neutral axis as the axis of no stretch and determing the projection 
of this axis onto the X-axis as the blade bends. This projection, then, is 
the spanvise or X location of the neutral axis in blade coordinates. The rate 
of change of this projection is the spanvise relative velocity and the second 
rate of change is the spanvise relative acceleration of the blade element 
neutral axis location or point. The motions are then transformed to the cen- 
ter of gravity to obtain the spanvise motion of the origin of the blade ele- 
ment reference axis. 


In Figure 27, the deflected neutral axis is shovs as a function of blade radius. 
The (i-1) and ith^ station are shovn. It can be seen from this figure that as 
X NA (i-l) approaches X^{ i ) , then the delta length of the blade 
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Rearranging this equation and summing from the blade root to the kth blade 
station yields: 
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NEUTRAL AXIS 


^BLn 



Figure 27. - Neutral axis vs blade radius. 

Likewise, 



( 11 * 2 ) 


S N ^(i) is simply the blade length to the ith station measured along the neutral 

axis and Y (i) and Z„ (i) are the Y and Z locations of the neutral axis in 
HBLn BLn 

the blade coordinate axis system for the nth blade. These displacements, 
along with their derivatives, will be defined later. First, however, by 
taking the first and second time derivative of X equation, the spanwise velo- 
cities and accelerations of the blade element neutral axis point are deter- 
mined and are given by the following two equations. 
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If Y rt _,(i) is the distance along the ith blade element chord line from the 
uha — 

blade element reference axis origin or center of gravity to the blade element 
neutral axis* then the blade element neutral axis motions can be written in 
terms of the blade element motions as: 
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Referring to Section k.5, the time derivative of the above equation is: 
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and likewise, the second time derivative is: 
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These three equations, then, define the Y and Z displacements, velocities, 
and accelerations of the neutral axis point used in the X equations and time 
derivatives. Also, the increments of spanvise motions due to the offset 
between the center of gravity and neutral axis are defined by these same 
three equations. This increment represents the motion of the neutral axis 
relative to the blade reference axis origin, therefore, the span motion at 
the center of gravity is determined by subtracting AX_,(i) from the spanvise 
motion of the neutral axis, or: "nLn 
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These equations, then, along vith the previous expressions for X and 2, de- 
fine the blade element relative displacement, velocity, and acceleration vec- 
tors required for the total inertial vectors which follow. 

1 *. 5.5.11 Blade motion in absolute coordinates. - To this point the blade 
element motion has been defined In terms of the blade axis or relative coordi- 
nates. Hie elements defined are: 


blade element relative displacements 
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blade element relative velocities 
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and blade element relative accelerations 
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Using the method of Section U.U.l, expressions in freestream (absolute) coordi- 
nates can be written for use in the equations of motion. The blade element 
velocity becomes: 
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The blade element accelerations are: 
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and matching rotation terms are defined in Section 4.5.4 in terms of rotor 
axis terms which are in turn related to the principal (hub) reference axis. 
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1*.5.6 Svashplate motion. - As shown in Figure 10, the svashplate reference 

axis system is defined with the Z-axis down. The motion of the swashplate 

reference system is defined by three generalized coordinate displacements, 

Z__, *, D , and 6™, which move relative to the hub axis system. 

SP OP or 

The rotations * sp and 0 sp are taken in the same order as shown in 
Figure 18 and therefore, from Section U.U.3 the angular velocities are: 
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where * sp 


is the rotational speed of the 


swashplate, and 


* 


SP 



(158) 


where is the rotational speed of the rotor. Note no coupling is pro- 
vided for shaft motion, the assumption being that swashplate motions relative 
to the hub due to shaft motions have been designed out of the system. 

As indicated before, the chosen svashplate axes do not rotate at the rotational 
speed * gp . Hovever, the total angular velocities reflect the rotational 
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rate Therefore, the total angular rates of the svashplate 

SF 

in svashplate axes are obtained with ♦ = 0. This gives 



(159) 


Bonrotating svashplate angular velocities, subscripted SP, are obtained 
by deleting *gp above. 


The svashplate angular accelerations can be similarly determined by 
evaluating the general expression at ♦ = $ = 0. This yields: 

Ox 



(160) 
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The vertical velocities and accelerations of the swashplate are simply 
defined as: 


2 0SP = 2 SP + 2 0H 


(161) 


and 


2 0SP = 2 SP + 2 0H 


(162) 


It is noted in these equations that the Z-axis motion is assumed to 
remain parallel to the hub Z axis. 

The swashplate angular displacements are obtained by integrating the 
angular velocities, or: 


■ A 


'SP - y * S P dt . ♦ s. 


( 163 ) 


and 


0 


sp / ®sp dt + e 

0 9 t=o, SP 


(16U) 


Likewise, the vertical displacement of the swashplate relative to the 
hub is: 


J 3P 


" / 2 


2 SP dt + Z 


t=0, SP 


(165) 
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U.5.7 Blade feathering notion . - The feathering occurring at the feather 
bearings. Figure 15, is taken to be the sun of the motions of the following 
dynamic and kinematic elements: 

• Swashplate - collective command 

• Swashplate - cyclic command 

• Blade bending to feathering couplings 

• Elastic pitch horn and associated components 
The total feathering response is: 

♦fd * 9 0 * *1S cos ^ 9 BLn * V - B 1S sin( *BLn * V * 


3 *Fn 

3A 1 A ln 


3*- 3#-. 

— £2. a + — — A 
3A 2 A 2n 3A 3 a 


+ — g 

3n 3Bp H °PHn 


(166) 


Velocities and accelerations are formed by differentiation, 
relations are: 

♦f- * % - *1S “ s( *BL„ * *R> - “is Sin( *BL„ * V 
♦ [a is slut*** ♦ *„) - B 1S cosft^ * * B )j 

. ilZB* „ 2*Sl i „ l . IflL X 

3A X A ln 3A 2 A 2n 3A 3 A 3n ae pH 6 PHn 


The desired 


(167) 
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and for accelerations: 



9 0 - ijs cosily * ♦„> - B 1S siB (* BLn ♦ ♦„> 


* 2 [*1S sl ” ( *BL» * V * hs coa( *BLn + Vl* 


*[* 


c ° s( *BLn 


+ B 1S 

sinC^ 

cosl *BLn 

* V 

- 3 1S 

sinU*. 

•L A + 

S *Fn 

A + 

3* f 

A 2n 

3A 3 

3n 

36 ph 


>]*, 

ili 


PHn 


a *Fn .. 

3A A ln 


(168) 


The coneaanded cyclic blade angles are: 



i 

I 



Si "*PH 

cos *PH 


*SP 

\e/ 

COStf’pjj 

- Si "*PH 


0 SP 


(169) 


where the angle is the pitch horn-swashplate connection lead to 

Jrii 

feather axis. See Figure 28. This angle is computed as a static value. 
It should be noted that some hub configurations carry the pitch horn 
toward the blade trailing edge. These configurations are entered in 
REXOR II by forming the supplement of 


i|i pH = 180 - il'pn (degrees) 


(170) 


This angle gives the correct modeling of the sense of rotation 
reversal with the trailing pitch horn geometry. 

The velocities and accelerations of the command cyclic are obtained by 
differentiation. 
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Figure 28. - Pitch horn blade feathering phase angle 


and 



( 171 ) 


( 172 ) 


The overall coupling (swashplate to feathering) gear ratio, d/e, is 
expressed as a static term plus a first-order collective correction. 



The collective is: 


6 0 = “ z sp /e 


(173) 


(17M 


The swashplate vertical motion, Z sp , is developed in Section 4.5.6. The 

value e is the static effective crank (pitch horn) arm about the blade 
feather axis. This crank length is entered as a negative number for a 
trailing pitch horn geometry to give the proper sense of collective for 
swashplate vertical translation. 
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Taking time derivatives: 


®0 * "^SP /e 


(175) 


and 


;; o - -V* (176) 

3 *Fn 3 ^Fn 

The blade bending to feathering coupling factors are — ? , tt — , and 

3A 1 3A 2 


3< ^Pn 

3A- 


for the first, second, and third blade modes. The blade bending modes 


are described without a torsion component; this allows freedom in varying the 
blade sweep, droop. Jog, or other geometric parameters without new input 
data for the blade mode shape. The torsion either is calculated separately 
along the blade proper or as a blade root component by pitch horn bending. 

The coupling factors are intended to add a feathering component to the blade 
mode which would exist even with no torsion or feathering moments . As such, 
they are in effect the 6,, a 0 » etc., coupling usually described in the 

literature. These couplings are usually determined as a function of the 
distance from the flap or inplane mechanical or vertical hinge to a pitch 
horn projection. 
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5. EQUATIONS OF MOTION 


5.1 Introduction 


With the coordinate systems tod transformation between systems well in 
hand, the development can proceed to the equations of motion. The 
development yields a set of second-order differential equations with time 
varying coefficients. These equations are formulated using the energy 
approach in a form credited to Lagrange. The solution to the system of 
equations is in the time domain by numerical integration. The result is a 
time history of the displacements, velocities, accelerations, and loads of 
the components of the helicopter modeled in detail, tod the program 
treats each blade separately. 

In the f< Hewing development a Lagrangian approach to system modeling is 
applied to a set of ppint masses and then extended to discrete masses and 
inertias. The result is a set of generalized mass and force expressions. 

In REXOR (Reference U) these expr Hons are programmed directly, element by 
element. In REXOR II extensive us_ is made of matrix notation both in descrip- 
tion and programming. The transition to matrix notation is given at the end 
of section 5*3. 


5.2 Energy Approach to Development of Equations of Motion 

There are two basic approaches to developing the equations of motion for 
a physical system. These are: 

• Vector summation of forces 

• Energy approach 

Given an equal set of conditions, limitations, and assumptions , both pro- 
cedures should result in equivalent sets of equations. The difference is 
in the ease of arriving at a complete set of equations. Note that force is 
a vector, whereas energy is a scalar quantity. Therefore, in dealing in 
terms of energy, less information regarding direction needs to be handled. 
Also the systematic nature of the energy approach reduces the risk of 
error. As stated by Lagrange (Mecanique Analytique, 1788), "The methods 
which I present here do not require either constructions or reasonings of 
geometrical or mechanical nature, but only algebraic operations proceeding 
after a regular and uniform plan". 

The starting point of this development is Lagrange's equation. It may be 
derived by postulating Newton's second law, or from Ila; .on's principle. 
Lagrange's equation may be written in the following fern.: 



(177) 
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where 

T is kinetic energy 
q is a generalized coordinate 
B is dissipation function 
U is potential energy function 

^ is the generalized force, derived from the virtual work, 6k, and 
is defined by the equation 



(178) 


Equation 177 will now be developed into the form as applied in REXOR II. 

This form bears a close resemblance to a force balance equation, but is 
derived from energy considerations. For clarity, the development is first 
shown for a set of discrete mass particles, then, in the section that 
follows, is extended to the distributed elemental masses of the REXOR II 
modeling and to the iterative solution scheme used. 

In a conventional manner the equation is formulated in terms of generalized 
coordinates. These coordinates are a function of time, and completely 
define the system. They sure genersilly not directly identifiable as a 
physicsd qusuitity. 

The physical parameters or elemental coordinates sure defined to be func- 
tions of the generalized coordinates and in turn a function of time. 

Consider a system to be composed of particles whose physical coordinates 
are a function of n generalized coordinates. For the ith particle: 

X i = X i^ q l’ q 2* ‘ ' ‘ ' q n ; (179) 

y i = y i (l l’ q 2’ * * * ’ V ^ (l80) 

Z i = 7 'i (q l* q 2’ ‘ ‘ * * V ^ (l8l) 

Note: a Cartesian coordinate set is selected, and used in REXOR II. However, 

the argument is true for srn arbitrary coordinate set. 
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The functional relationship of the physical or constrained coordinates and 
generalized coordinates yields: 


6x. 

1 


3x. 3x. 

— -i 6q + r- — 6a ♦ 

3q i 1 ^ 2 


3x 


(182) 




6q x + 


h 2 




( 183 ) 



3z. 3z. 

— - 6a + — — 6a + 

*1 q l 3<L, *2 


3z. 


(18U) 


The time dependence is lap’ \cit in the inc rem ents of the generalized 
coordinates. The equation is strictly true for infinitesimal increments . 

In REXOR II the generalized coordinates are distinct frc*n physical coordinates 
in the main rotor blade descriptions. Here the generalized coordinates 
are blade modal variables. The modal variables represent tangible deflec- 
tions of the blade frees a reference position, and as such sue small but 
not infinitesimal variables. 


As the variables are a function of tine: 


3x. 3x. 

X i * 3q x q i + 3q 2 q 2 + 


3x. 


( 185 ) 


3y. 3y. 

y • - % ♦ ■^~ = ' 4-. ♦ 

l 3q x 1 aq, 2 


*i - 




( 186 ) 


3z. 3z . 

Z i ~ q i + q 2 + 


3z. 


( 187 ) 


In terms of the ith purticle the kinetic energy for the system may be 
identified as: 


is 

L I ( 2 .2 . 2\ 

2 "i \ X i + y i + Z i ) 


(188) 


i=l 
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Toward the particular formulation of Lagrange's equation used in REXOR II, 
the first two terns of the previously stated form. Equation 177, are 
developed: 



(189) 


Performing these operations for the ith particle case and the rth gen- 
eralized coordinate and summing over the system yields: 



A useful math operation of cancellation of the dots is developed prior to 
proceeding. Recall : 



3x. 3x. 

+ T — ~ <5q_ + 
3q l 1 3q 2 ^ 


Then also 



3x. 3x. 

X - 




3x. 

+ ^ 4n 


or 


3x. 3x. 

l l 



(191) 


( 192 ) 


(193) 


This is also true for y and z and for the double dot terms in x, y 
and z. 


An operation to reverse the order of spacial and temporal differentiation 
is required. To show this the time derivative of a partial is taken as 



(19b) 
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Next the special derivative of 


is given as 


3x. - fix. 3x 

— - = — — - 4 ♦ — - i 

*1, 3q r\ 3q l ^ q 2 92 


♦ - 


3x. \ 


Kov since 


x. 

1 


x i l V ’a- 



(195) 


(196) 


the order of spacial differentiation is reversible 


3q r 3q s 3q s 3q r 

and hence 

3x. 

l 

3q 

r 

Similar relations exist for y. and z. . 

l l 

Proceeding on with the kinetic energy terns: 



(197) 


(193) 




3.x. 3y . 

l - l 





(199) 


Using the relationship for cancelling dots in partials, reversing the order 
of differentiation and cancelling terms gives 



( 200 ) 
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Then from Equation 1T7, Lagrange's Equation in constrained coordinates 
vith point masses becomes 



(201) 


Also, in the same vein of defining the generalized coordinates, the 
relationship between the elemental and generalized forces can be developed. 
This relatic nchip is developed from the definition of virtual work on a 
particle as the scalar product of the applied force and an infinitesimal 
displacement . There "ore for the total system of N elements. 


6V 




6x. + F 6y. 
1 y. 



( 202 ) 


Using the definition of 


from Equation 178 gives: 



(203) 


Substituting Equation 203 into Equation 201 yields the final form of the 
Lagrange energy equation in constrained coordinates for point masses, which 
is in the form from which the REXCR II Equations of motion are developed. 
Making this substitution and rearranging the equation yields 





3B 3U_ 


0 


(204) 
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The above equation is the basis for the entire derivation of the equations 
of notion of REXOR II. note that this equation is written for discrete ele- 
ment masses and discrete forces. Also, at any instant in tine all of the 

ingredients required to define the elemental accelerations, x^, y^, z^, 

are not known. Specifically, the generalized coordinate displacements and 

velocities, q^ end q^, are kncvn at any instant in time but the generalized 

coordinate accelerations, q^, remain to be determined at the time the 

elemental accelerations are computed. 

The following section presents the manner in which the foregoing equation 
set is adapted to the REXOR II numerical solution to solve the equilibrium 
equations or equations of motions for the generalized coordinate accelera- 
tions. This development is first presented in the simpler form, for 
clarity, for discrete mass elements and forces and then in expanded form 
to include elemental distributed masses and applied moments. 


5.3 Itera tive Concept and Equation Set Solution Method 

Given a set of equations as developed in the previous section, the next 
step is to establish a method of solution. The solution process is 
defined as solving the equation set for the accelerations, integrating 
the accelerations for updated velocity, and position; then substituting 
the integrands back into the equations to determine new values of 
accelerations. 

•Hie first step of the process is to define exolicit.lv t.he 

from the equation set. In the process of implementing the REXOR II equations, 
it is desirable to handle the accelerations as an estimated plus a correc- 
tive term. In generalized coordinates then we can write 


r V 


'V 


’ V 

• 

► = < 

• 

► +■ 

• 

A 

NEW 

A 

CORR. 

• 

AJ 


(205) 


This operation proceeds on a sequential time basis. For each increment 
advance in time, the previous ’NEW' becomes the ’OLD'. In REXOR II, the time 
increment corresponds with a step azimuthal advance of the main rotor 
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blades. However, this need not be the case. The ’NEW' accelerations must 
be used in the numerical integration process to define the generalized 
coordinate velocities and displacements. But if some form of a predictor 
on accelerations is used then the 'OLD' would be this predicted value and 
in this case it would be an estimated, *EST', value. 


Using the notation 'OLD' and 'EST' interchangeably the linear elemental 
accelerations can be written at time t as 



(206) 


where the estimated accelerations are determined using the generalized 
displacements and velocities, q^ and q^, at time t, and the generalized 
coordinate accelerations q^, either estimated or from one previous time 
step in the numerical integration process. 

Then, at any given instant in time where the ’EST' elemental accelerations 
are thusly determined, the corrective elemental accelerations, (x, y, z CORR 
can be written as a function of the generalized coordinate corrective 
accelerat ions . 

Or 


r ' 

X. 

1 


r 3x. 3x. 

"5 — + . . . + ■— 1 d 

3q l X C0RR 3q n "CORR 

*i 

► = < 

ay i . 9y i .. 

— q + • * * + — — d 

. q l X C0RR 3q n ^CORR 

z. 

1 

CORR 

3z. 3z. 

- — q + • • • + - — d 

3q l X C0RR 3q r. ^CORR 

V . 
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or 



Nov mailing the substitution of Equations 203 and 207 into Equation 201 
from the previous section and rearranging terms yields the Lagrange 

equation for the q^ coordinate in terms of the estimated elemental 

accelerations and the corrective generalized coordinate accelerations. 
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The equations of motion for the system can now be combined and presented 
in matrix form. 




( 210 ) 


where the matrices, M^, C rfc and K rfc will be defined in the following 

discussion. However before proceeding with this, Fquation 210 is now 
rearranged into the form actually used in the numerical process in REXOR II . 
The equation is solved in terms of the corrective accelerations. 
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The correction terms come from an inversion (or simultaneous equation, 
Cholesky method) operation on the model equation set. 




'V 




'V 

C rk 

• 

y + 

K u 

rk 

< 

• 

DAMPING 



STIFFNESS 



MATRIX 



MATRIX 


* 




_ _ 




( 211 ) 


As indicated before estimated accelerations in physical coordinates come 
from the 'EST' or 'OLD' generalized coordinate accelerations and the cur- 
rent generalized coordinate velocities and displacements. The integration 
part of the solution operation supplies the (q) and (q) data. 


q = 



dt 



( 212 ) 


The whole package operates in a cyclical fashion, as shown in Figure 29- 
Arranging the solution sequence as such gives it some important attributes 
and advantages. 


First, to determine the corrective acceleration, the inverted mass matrix 
premultiplies the difference of applied and estimated reactive forces 
represented by the quantity in the large brackets on the right-hand side 
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NEW ACCELERATION 
f> REPLACES OLD 
ACCELERATION 


FORM NEW AlCEI E RATION 
FROM OLD AND CORRECTION 
TERMS 




INTEGRATE NEW 
ACCELERATION TO 
FORM VELOCITY AND 
POSITION TERMS 


ADVANCE 

TIME 



V 

SUBSTITUTE ACCELERATION. 
VELOCITY AND POSITION 
DATA BACK INTO EQUATION 
SET 


SOLVE FOR 
ACCELERATION 
CORRECTION TERMS 



Figure 29. - Equation solution loop. 
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of Equation 211. With the usual, small, integration steps these 
differences will be relatively small. Therefore, inaccuracies in the 
mass matrix or its inversion process only slightly affect the total 
acceleration determination. This means approximations and simplifications 
to the mass matrix are acceptable. In some instances, a diagonal mass 
matrix will give convergence to the required solution. 

Second, as will be shown in the Section 5.^, (blade equations section), 
carrying the running acceleration in elemental coordinates allows for the 
simple separation of the centrifugal and structural stiffness of the rotor 
blades which has important advantages which have been discussed. Also, 
the aerodynamic loading terms, already by nature in physical coordinates, 
are easily accounted for. 

In the actual application of Equation 211 to REXOR II, distributed elemental 
rigid body masses are associated with each coordinate point and applied 
moments in addition to forces at each coordinate point are accounted for. 

Referring back to Equation 203 the generalized force, Qy, from virtual 

work can be simply written in the following form to account for applied 
moments at each of the ith grid points as 



The terms of Equation 200 in Equation 20U can be developed for the dis- 
tributed masses by going back to the elemental acceleration equation. 
Equation 13 of Section h.l+.i which is repeated here, in a rearranged 
form, for clarity of this development. 
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For distributed masses of a rigid body with coordinate point and system 
embedded in the body: 


x = y = z = 0 


1 215) 


and Equation 2lU oecomes: 


! -x(r 2 +q 2 ) + ypq + zpr| | zq - yr 
xpq - y{r 2 +p 2 ) + zqr >+ <-zp + xr 
xpr + yqr - z(p 2 +q 2 ) I I yp - xq 


(216) 


Now, remembering that for a point mass. 


and 


3x _ 3x 

8< i r ~ 3 5 r 


ay = 3jL 

3 9j. 


3z _ 3z 


3$ _ 3p 

3< l r " \ 


36 _ 3q_ 


3<fc 3r 



(217) 

(218) 

(219) 

( 220 ) 
( 221 ) 

(222) 
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The total partial derivatives relating the motion of the coordinate point 
end set imbedded within each elemental body and the motion of the 
generalized coordinate becomes 


3x _ 3x o ae 

8q r 3q r ' y 3q r 2 3q r 


(223) 


ay _ 3y o t 3<j> 

x 3^ “ 2 3^ 


(22U) 


3z _ ^0 _ 38_ 3$_ 

9q r 3q r ~ X 3q r ^ 3q r 


(225) 


where on the right side of these equations, x, y, and z represent the 
location of the distributed masses within the rigid body elemental mass, 

and Xq, y^ and z^ represent the motion of the mass element reference 

point. 

For each Jth_ coordinate of the system, the elements of 

Equation 200 can be written by substitution of Equations 21 6 , 223, 

22b and 225 • This gives 
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Expanding and identifying mass moment and moment of inertia terms: 
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and finally collecting and grouping terns yields the final and complete 
definition of the tens of Equation 2xi x'or the estimated elemental 
accelerations . 
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where in this case the summation 


represents suimnation over the jth 


rigid bcd;r element. With this in mini, substituting Equations 228 and 
2_3 back into Equation 201 yields the complete form of the Lagrange 
energy equation in constrained coordinates with distributed elemental 
masses and forces from which the REXOFt II equations of motion are all devel- 
oped. Also including these terms as veil as the moment terms of Equa- 
tion 213 in Equation 211 yields the final form of the equation as used 
in REXOR II. This form will be presented following the development of the 
generalized mass damping and stiffness matrices. 

From Equation 209 it is easily seen, by examining the coefficients of 
the corrective accelerations that generalized mass matrix elements, 

M , , can be written as 


hx. 3x. 3y. 3y. 3z. 3z.\ 

rlt " 3q k 3q r 3q k 3q r 3q k/ 


( 229 ) 


This equation is for point masses. Actually, as discussed earlier, the 
REXOR II equations model a set of distributed masses characterized by an 

overall mass, center of gravity, and moment of inertia values. As shown 
in the previous section, extension to the distributed mass form is made 
by describing the particle absolute coordinates in terms of the position 
of a relative coordinate set 'n inertial space and the particle position in 
terms of this relative set as developed in Section U.U. For a rigid body 
the associated relative set and the particle associated with the body main- 
tain a fixed relationship. The summing over the particles of the system 
then becomes a sum over products of masses and lengths yielding mass 
moment and moment of inertia terms. 

The mass elements can be a ./eloped by substituting the partial derivatives 
developed in the preceding discussion. These partials describe both the 
motion of the mass element reference and also the distributed masses within 
the rigid body elemental masses . 
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Substituting these partials, Equations 223, 22U and 225 in the 
generalized mass expression. Equation 229 yields: 
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and using moment of inertia and mass moment definitions: 
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and is identified as a generalized mass. For orthogonal systems M is 

r k 

zero except for r = k. The development of REXOR II is mostly nono rthogonal 
coordinates, therefore, the generalized mass matrix has many off-diagonal 
terms . 

Similarly, terms can be developed for the strain (potential) energy and 
damping functions. 
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Similarly for damping: 
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The stiffhess and damping matrix terms In REXOR II are defined with reference 
to relative coordinates; which parallels the physical configuration. The 
coordinates used with these terms then should be on a relative basis. 

This statement at first appears to be contradictory to the premise of the 
equation development . However, if these matrix terms were defined on an 
absolute basis the terms other than those associated with a relative 
motion would be identically zero. The integration of the accelerations 
produces changes in velocity and position. These changes with the proper 
starting reference are the relative coordinates and velocities. 

Equation 211 is now repeated here in a slightly expanded form to include 
the effect of applied elemental moments. Equation 213, and distributed 
elemental masses. Equation 228. 






(236) 
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Even though vas defined as the generalized forces of the system, for 

the purpose of further development and of the application of the above 
equation in the KEXCR II analysis, each line in the large brackets on the 
right side of Equation 236 will hereafter be referred to as a generalized 
force or a generalized delta force and will be referred to by the symbol, 
FR, in the following development. 


In REXOR II use is made of matrix notation to produce compact partial 
derivative, generalized mass and generalized force expressions. The 
partial derivative set 

ax ay 3z ae at 

a^’ aq^ aq r ’ aq r * aq r ’ aq r 

is replaced by 


3t 

3< *r ' 


Usually the generalized coordinates q are grouped as three linear plus 
three rotational motions. Hie full partial derivative for system coordi- 
nate A and generalized coordinate set B is: 



The generalized mass expression (Equation 231 ) can be rewritten as: 
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The general! ied delta force, Ao , can be expressed In matrix notation for 
contribution of coordinate set A as: 
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5. ^ Overview of Rotor-Blade Model 

Many elements of the rotorcraft can be directly modeled following the 
methods developed in Sections 5 through 5*2 and systematized in Section 5*5* 
However, there are enough special considerations and concepts involved in 
modeling the individual blades and combined rotor to Justify a separate sec- 
tion to address these topics. 

5.U.1 Concept of modes . - The basic textbook principles governing solutions 
for eigenvalues (natural frequencies) and eigenvectors (mode shapes) for 
systems of several degrees of freedom can be applied to those of many degrees 
of freedom. For each independent degree of freedom there is an additional 
natural frequency and mode shape. 

Free vibrations of continuous systems such as beams, or for example the 
helicopter fuselage, or rotor blades, are generally analyzed mathematically 
by reducing the system to a system of discrete masses and elastic 
constraints . 

5.U.2 Blade bending - modal variable . - The blade is a twisted rotating 
beam and its analysis requires considering the coupled flapwise-chordwise- 
torsional response of the blade. For the REXOR II analysis, coupled flapwise- 
chordwise mode shapes are used, upon which is superimposed one of a number 
of torsional response representations of varying complexity (Sections h.3.h t 
5.U.T, 5*6.5, and 5*6.6). 

If one applies generalized coordinates, each blade mode in the analysis 
may be treated as a single degree of freedom. The generalized coordinates 
are called normal coordinates for the special case when the modes are 
orthogonal, in which case the generalized mass matrix reduces to a diagonal 
matrix, as does the generalized stiffness matrix. 

The REXOR II analysis uses blade modes calculated for the blade at a fixed 
rpm, fixed collective, j «d ip an unsvept, unconed orientation. Since the 
program allows for variation of all of these parameters, which is accounted 
for in the overall REXOR II analysis, the predetermined modes become non- 
orthogonal as used in the program. Thus, blade motion is effectively 
described by a set of modal variables, each representing a characteristic 
frequency, and a set of modal coefficients that describe the relative 
amplitude of oscillation for each blade segment and each frequency. 

Since the modes are nonorthogonal , we will find in REXOR II, as would be 
expected in such a case, off-diagonal coupling terms in both the gen- 
eralized mass and stiffness matrices. It can readily be shown in cases 
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where generalized or normal coordinates are applied, that relatively few 
modes need be taken to define accurately the time-history of blade deflec- 
tion. This assumes that the primary frequencies of excitation fall within 
the range of modes considered. 

5.4.3 Blade mode generation . - The blade modes can be determined by any 
appropriate classical method of analysis for coupled flapwise-inplane bending 
beams. The only requirements is a cantilever (hinge or hingeless) boundary 
condition for the modes and that the terms included in the homogenous part 
of equations 28 and 29 of Reference 5 be accounted for. These equations are 
repeated here for convenience. Flapvise: 


£(EI l cos 2 B + EI 2 sin 2 6)v" + (EI 2 - EL^ sinS cosB v"J" 


- (Tv')’ - ft mv + mv = 0 


(239) 


and inplanevise: 


£(EI 2 - EI^) sinB cosB w" ♦ (EI^ sin 2 B + Elg cos 2 3)v"J l 


- (Tv' ) ' - ft mv + mv = 0 


( 21 * 0 ) 


5.4.4 Modal coefficients . - Several additional points need to be made 
regarding modes in order that the equation development be properly understood. 
First, the same modal coefficients apply to the first and second time 
derivatives of the function, since 



f ( x ) 


(241) 
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Then 
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Second, the motion is not necessarily confined to one direction. A given 
modal frequency may excite or couple with motions in other directions. 

For example: 
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5.4.5 Independent blades . - In REXOR II the blade motions are computed and 
tracked individually. 5ne set of equations operates on a blade in BLE 
coordinates as explained - Section 5.4.11. The result for a time step is 
stored in BLn coordinates for that blade. The operating set in BLE 
coordinates then performs the computation; for the next blade in turn. 


5*4.6 Blade element aerodynamic forces - overview . - The functions F. 
Fy^, F^, momen t terms from Section 5.2 are primarily aerodynamic 


*i’ 


loads for the blade equations. These loads are derived from blade inertial 
velocity (equivalent to air velocity) and table lookup aerodynamic 
coefficients as given in Section 6. 
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5 .4.7 Blade torsional response. 


5. 4. 7.1 Pitch horn bending . - Several alternate approaches to modeling blade 
feathering dynamics exist in REXOR II. One approach is to assume the blade 
is torsionally rigid, and that the flexibility is in the pitch horn. 

5.U.7.2 Quasi-static blade torsion . - The blade pitch horn bending descrip- 
tion is improved by the addition of a blade twist dependent on the moment 
loading. This quasi-static torsion is computed by integrating the blade 
pitching moment times the torsional flexibility from tip into the root. 
(Developed in Sections 4.3.4 and 5.6.5.) 

5.4.7. 3 Dynamic blade torsion . - A third approach to blade torsional response 
in REXOR is an uncoupled torsional mode which operates as additional blade 
twist. This material is developed in Section 5.6. 

5.4.8 Radial integration . - For each element of a rotor blade the equations 
of motion are formed per Section >.2.9. As briefly touched on in Sec- 
tion 5.1*. 6 these data are formed in L>LE axis. These elements are summed to 
total equations for each blade in BLn coordinate at the blade root. This is 
explained in Sections 5.6.3 and 5-6.4. These blade root simulations are also 
used in the fuselage axis (Section 5.8). 


5.5 Equation System Development 

5.5*1 Reference to base operation matrix . - The equation of motion, as 
developed ir. Sections 5.2 and 5-3 and as presented in most general form by 
Equation 236, may be given in abbreviated form as 

{ * r ] CORR " [" rk ] KJeST (2l,7) 

The M r jj' s are the generalized mass matrix elements, the F r 's are the 
generalized forces, and the q r 's the generalized coordinates or degrees 
of freedom. As explained previously, the F^'s are the complete set of 
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external forces and internal reactions computed with estimated values of 
the accelerations, <i rEST ’ s » at the next time point. The Si-corr ' 8 are 

then corrections to the estimated values. 

The generalized mass, M ^ is developed in Section 5.3. The generalized 
force may be expanded as (using the point mass form): 



+ F + F. + F (2U8) 

FR A c 


The inertia, damping, and elastic terms are developed further in Sec- 
tion 5.3 (see Equation 236). The friction force F pR , the aerodynamic 

r 

external forces F. , and the pilot control forces F are described as 
A c 

r r 

needed. Note that the potential energy and dissipation terms have been 

directly included in the force expression. Where the stiffness and damping 

matrices are simple diagonals, this is done. In the case of the blade 

equations the distinct stiffness and damping matrix form (Section 5.3) is 

computed before combining all the applied forces, internal reactions, 

stiffness and damping terms into an overall force. 

5.5.2 Organization by degrees of freedom . - In developing the equations of 
motion there are three types of ingredients needed: 

• Generalized masses 

• Generalized forces 

• Partial derivatives (used in both of the above items) 

The equation development ceuj then proceed witn these ingredients along one 
of two lines of organization. 

• For every major rotorcraft piece (fuselage, rotor, etc.), compute 
all the ingredients and sort according to degree 01 freedom for 
equation use. 



• For every degree-of-freedom group, sort through the rotorcraft 
pieces for applicable ingredients. Sorting is minimal because of 
close association of degrees of freedom and component parts. 

The latter development is used here. The degrees of freedom modeled are 
given in Figure 11. 

The following subsections will describe the appropriate partial differen- 
tiations, the generalized masses, and the generalized forces in detail. 

Each generalized mass couples the inertia of one generalized coordinate 
with another or itself. The algebraic equations for each generalized mass 
will be given only once. If the reader cannot find a particular mass 
element under one subsection, he should look into the other subsection 
relating to the coupled generalized coordinate. 

5,5.3 Partial uerivatives . - Tne generalized masses and forces use partial 
derivatives which cescr oe the variational motion of each physical mass element 
in rectangular coordinates relative to the motion of each generalized coordi- 
nate. The partial derivatives required are determined from the generalized 
mass and force expressions for distributed masses of Section 5*2. The partial 
derivatives are easily constructed from the coordinate transformations which 
have been developed. 

In developing the motions of a physical mass element relative to a gen- 
eralized coordinate, a number of transforms may be used. These can be 
categorized as either linear or Euler axes transforms which either displace 
without rotation or rotate without displacement. The overall partial will 
be the product of partials associated wfth each of these transforms. The 
typical form of these partials will now be illustrated. 

To obtain the partials, the equations relating the velocities are obtained 
first. Reviewing Sections k.4.1 ana h.h.2, the velocity relations of 
interest are restated. For linear transforms: 
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and for Euler transforms 






The parti&ls of interest are conveniently organized into 3 by 
They are for the linear transform: 



► 

X 



X 



X 

' - 



a 

v 

t 

3 

V 


3 

v 



M 

3X n 

0a 

I 


3Y_ 

0a 

I 


0a 

I 

► 


L J J 

- 

z 

b 


z 

k. 

. b 


z 

b_ 



- 

X 



X 

l 1 

’ 

X 

-i 



a 

3X ‘ 
a 

Y 


3 

3Y < 
a 

Y 


1st- 

a 

Y 


= 

H 

- 

Z 

k. 

b 


z 

k a 

b 


Z 

• . 

b 



m 

X 

l f 

> 

X 



X 

- 


0 

3 

3*a - 

Y 

► 

4- 

Y 

► 

3 

Y 

> 

= 

-Z 


Z 

4 

b 


Z 

d 

b 


z 

d 

b 


Y 


10 0 
0 10 
0 0 1 



litU 


(239) 


► V (2U0) 


matrices . 


(2U1) 


(2k2) 


(?M) 




H»5 




► « 

$ 

1 

a 

r 

$ 

1 

■■ ' 
* 



cos^cosQ 

siniji 

o' 

a 

3 V 

e 

' 1 k 

1 b 

e 

1 !*■ 
I b 

e 

> 

= 

-sin^cosd 

cos^ 

0 


$ 

fc 4 

* 1 

* 

b 4 

» 1 

* 

» 4 

b 


0 

0 

1. 


( 250 ) 


For Euler angles defined in reversed order or reversed sign the last matrix 
will differ. Note by inspection that rotary to linear derivatives such 

36 

as — are all zero. The derivatives can be strung together to get motion 

in a third axis c relative to motion in axis a. Abbreviating the 
matrices: 
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assuming in geneiol 



The abbreviations used are 


r = {X, Y, Z} (257) 

and 

;={♦, e,*> (256) 


Transformations may involve linear and rotational operations. Partial 
derivatives shoving the combined operations may be generated using a 
linear transform from set a to b, followed by a rotation from b to c. 
This sequence gives: 

To 
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T— 1 = 0 

U r aJ 

BS1-M 


5-5-4 Generalized masses . - As discussed before, the helicopter is assumed 
to consist of a finite number of mass elements. They are the 

e fuselage 

• tail rotor 

• engine rotor 

• svashplate 

e fixed hub {all parts that do not feather) 

• k mass elements on each of b blades. 
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The reader should realize the mass matrix is syanetric from the definition 
of Equation (209) and interchange of the order of differentiation. 


M. * H 

xr 


rk 


(263) 


Only the elements in the diagonal and the upper right triangle will be 
given in the following sections. 

Each of these mass elements must be summed for each of the generalized mass 
matrix elements. Each mass is handled with the distributed mass relation 

of Section 5-2. Fortunately, only the fuselage requires the full equation. 

The center-of-gravity terms drop out if the mass motion is determined at the 
center of gravity. This situation is true for the blade line which passes 
through the center of gravity of the blade section mass elements . Only the 
fuselage, the svashplate and the shaft/transmission have reference axis origin 
off the center of gravity. Another simplification is that cross products of 
inertia exist only for the fuselage and the shaft/transmission. Each blade 
mass element is considered to be in the shape of a rod lying along the chord 
at the blade station in question. 

Certain small terms and factors are dropped from the generalized masses. 

As discussed in Section 5.3, the equations of motion are solved for small 
incremental corrections to the accelerations. With this formulation the 
masses car. tolerate approximations as contrasted to the generalized forces. 

5-5.5 Generalized forces . - The equation formulation. Section 5.3, requires 
that precision be used in compiling the generalized forces per Equation 247, 
expanded per Section 5.3, Equation 236, to include rigid body distributed 
mass elenents. This formulation includes for each degree of freedom: 

• Suianation over all mass elements of the mass times inertial 
acceleration times paitial derivative. (Section 5.3 expression 
for distributed masses.) 

e External (aerodynamic) loadings times a partial derivative. 

e Potential energy and damping terms or assembled stiffness and 
damping terms with partial derivatives (Section 5.3). 

For some degrees of freedom the applicable mass elements and the total 
integration are directly written as final results which can be verified 
by inspection. Degrees of freedom that properly include summation over the 
main rotor blades involve some extensive numerical integrations and com- 
plicated coordinate transformations. 
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5 .6 Blade Bending and Torsion Equations 


5-6.1 Blade radial summation . - The contribution from all the individual 
blade sections are summed to give the blade generalized masses and forces. 
These are given for blade root, bending, feathering, and torsion motions. 

The blade root values are then transformed to the final degree of freedom 
variables by partial derivatives. The summation is carried out over all 
elements of the rotorcraft, including the independent blades. Due to the 
relative isolation of one blade's modes from another, only the 4 by 4 submass 
matrices along the diagonal of the 4b by 4b rotor matrix are filled, where 
b is the number of blades. 

5-6.2 Partial derivatives . - The generalized masses and forces utilize 
partials relating the X, ¥, Z, $, 6, and linear and rotary notion of each 
blade element to the blade bending, blade torsion, body, rotor, and swash- 
plate degrees of freedom. Only the blade bending, torsion, and feathering 
partials are derived in this section; the blade partials for other degrees 
of freedom are to be found in their respective sections. 

As developed in Section 4.3.4, the blade torsion may be modeled either as 
a pitch horn bending or an uncoupled dynamic torsion mode. For the former 

a *Fn 

case the partial — — is a blade spanvise constant multiplier to sussaa- 

PHn 


tions which couple in the feather angle. 


In the latter case. 
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as, 
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is a 


function of span and blade number. 

The first partials to be considered are those relating motions at any 
point i or. the blade to the rigid body notion of the blade root. These 
partials are: 
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8ote that 


X BLE’ Y : 


and are expressed in blade root coordinates; 


while and <*gjjg are in tertas of blade element axes aligned 


with the blade element principal axes. 


Next consider the blade Y and 2 bending response with respect to the 
blade bending modes. A number of equations can be used to develop the 
required expressions. The velocity equations from Section 1.5*5 are 
selected for ease of analysis. Using cancellation of the dots 
(Section 5.2): 
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where 
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The angular derivatives with respect to the blade bending modes are also 
constructed in the velocity form. 
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Note that the angular derivatives, being applied to local segments, are 
presented in BLE axis. Referring to Equation 132: 
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Also note that in the same context and argument of Section U.5-5 the 
feathering axis slopes, Y'^ and Z have been neglected in the above 

angular partials. 

Derivatives with respect to the blade feathering are also constructed 
using cancellation of the dots. 
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Similarly, for angular motion: 
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The partials developed with respect to are used directly in swash- 

plate and rotor summations as well as some of the following mass and force 


terms . 


Some terms require a further compounding derivative. 
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the case of the pitch horn bending torsion option. Talcing the deg ■* of 
freedom to be pitch horn angular deflection about the feather axis, the 
constant is approximately 1/e. For the dynamic torsion option a different 
set of mass formulations is used in terms of BLE axis, obviating the 
need for the compounded derivative. As indicated in Section 4.5.5, 
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Blade X motions must now be accounted for. The assumption of the neutral 
axis as the axis of no stretch is discussed in Section U.5-5 and the 
derivation of the X motions shown. The equation for the partials in 
BLn axis for a point on the neutral axis is taken from the formulation for 
the X velocities: 



The program data, however, is at the blade center-of-gravity axis. The 
transfer is : 
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and 




A simple partial derivative is also needed when is defined as dynamic 

torsion. Since torsion occurs along the bent and twisted blade lire, in 
blade element axis BLE, only the vertical cr normal to chord motion of the 
shear center is of interest, hence. 
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38 PHn 


is program input for the torsion node shape. 


Use is nade of a blade node to blade feathering partial derivative array to 
produce a compact development . This array is simply by definition: 


ar BL 

d*Fn 


[°] 


0 0 0 ' 1 


(287) 


Where | T Bl| = | A ln* A 2n’ A 3n* *Fn|* 

Additional partitl sets which are used to expedite the b&ss and force 
expression development are: 
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$.6.3 Generalized masses . - The blade generalized masses in conjunction 
vith partial derivatives couple blade feathering, blade torsion, blade bend- 
ing, fuselage motions. All the blade generalized masses that assumed to 

3X BLE 

exist are given in the following table. As mentioned before, 3 “ is 

Tn 


assumed zero and 


3 * 


BLE 


H- 


Fn 


is assumed one in the program. 


although it is given 


in the table, 
axis I—y. 

BLE 


The blade has a rotary inertia about the center of gravity 
The blade also has inertia I__ about a vertical axis. 


Table 1 lists all the terms coupling rotary motion at the blade root, 

M and similar terms. However, not all listed are used, as certain 

*BLn eLn 

approximations are made in developing the principal axis generalized masses 
which reduce the number of blade coupling generalized masses needed. Since 
the mass matrix operates on the acceleration error term rather than the total 
acceleration, these approximations do not detract frcr the validity of the 
results produced. 
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TABLE 1. - Continued 
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To save computation time some of the masses generated by the blade inter- 
gration and summation process in Table 1 are saved as a pseudo mass associated 
with a fictitious rotor coordinate (R). Thus for any given time step the sum 
of the rotor blades can be treated as an equivalent mass matrix: 



The couplings of the pseudo rotor coordinate to blade and hub coordinates are 
also formed as an intermediate step to save repetitious blade integration. 
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Where the array 


is the feather and bending coupling terns for 


M. A *X A A 14 *.* 
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In addition to the rigid body motion blade coupling matrix, M_ , the 

oL-BLj 

blade mode coupling matrix, W BL _ BL » wil l be used in the subsequent 


development 
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5.6. 1 * Generalized forces . - The development herein proceeds by first deriving 
the equations for the loads on an individual blade element. The blade element 
loads are composed of aerodynamic and inertial components conveniently found 
in either the blade root axes BLn or the blade element axes BLE. The loads 
will be summed in Bln axes with the appropriate transformation. The desired 
equations are: 
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The aerodynamic loads are in BLE axes alignment about the blade reference 
datum line which is the quarter chord. A transfer through the distance 
Y__ is made to the aerodynamic moment. To put the data on a common 
BLE 



basis with dynamic terms. The blade aerodynamics is detailed in 
Section 6. Since only the blade section pitching moment is considered, 

XlE ' "“BLE ‘ °‘ Note the blade element is assumed configured as a 
chordvise rod for inertia; hence 


( I * I and I s 01 

Xle zz ble ^ble ) 
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A number of blade summations are desired. All will be made in BLn axes 
along the center-of-gravity axis. The loads at the principal reference 
axes and for rotor tilt make use of the blade root shears and moment. 
These are simply the sum of the k total blade elements. 
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and likewise for root moments. 
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The simulations illustrated above are for the total inertial and aerodynamic 
components. In a similar manner, the blade root aerodynamic loads are 
derived. The blade root loads are summed over all the blade to give main 
aerodynamic loads for downwash computations (Section 6.2.2) in the manner 
the total main rotor loads are found in Section 5.9. 
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Total main rotor root loads are formed from the blade root shears and moments 
Using the pseudo rotor coordinate: 
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Feathering moments are used by the svashplate equations of motion. These 
moments are: 
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where 
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Only the X component is used by the program. Thf equations above trans- 
fer the sunned blade loads to the inboard bearing, hen transform them to 
feathering axes. Using the blade root loads is • i ect when one recalls 
that the blade is defined as those portions that ■. « feathered; the fixed 
hub is excluded. 

The blade bendiig generalized forces are now presented. 1 y are: 
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The potential energy is given as : 
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ROOT 

are inputs calculated external to the program from a be-^ing beam model. 
El^y and EI ZZ are the flapping and chord stiffness abc t axes aligned 


with the blade element principal axes. The chord and flapping moments, 

((v . and M reflect the contribution of the bending moment from the i 
T i Z i 

(or j) mode. The integration goes from root tc tip. The K’s are 
evaluated for whatever normalized modes are used as program input. 


The last equation can be derived from the Bernoulli-Euler law for bending 
beams : 
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where r is a radius of curvature. 


The strain cergy is 


'/root 2 (*» 4S * ^ «s) 
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Substituting in from the Berne uili-Euler lav and noting that dS * r v d9 = 


r„d$. 



(330) 


Partial differentiation gives 
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Considering the moments as a linear sum of components fro® each bending 
mcde, one has: 
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and likewise for Si,. Then, by substitution, the desired equation is obtained. 

I * 

5. 6. 1.1 Blade motion dampin g. - The blade notion damping is modeled by struc- 
tural damping, aercdynar.it camping, and a m^hanical aampt-r for lead- lag 
motions. The aerodynamic damping is accounted for in the aerodynamic blade 
loads developed ir lection 6.2. 


The structural factor is assumed proportional to the spring ’-ate components 
of "cuaticn 327. The coefficients K , can be directly identified with the 

coefficients " developed in Section 5-3, Equations 231 and 236 The required 


camping coefficients, C^, are then tfc.- product of the proportionality factor, 
T.., and '< The structural dancing component is then; 
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A lead-lag mechanical damper is usually a rotary or linear motion device em- 
ploying elastomeric, dry friction, or viscous energy absorbing mechanism. 

The device is oupled about the blade lead-lag hinge or point of lead-lag 
motion slope by a linkage array. REXOR II models a rotary viscous damper 
mounted about a given blade location for which the inplane slope is specified 
as a function of the blade modal variables. Pbr an articulated blade, this 
will simply be the motion of the lead-lag hinge. The damper has a pressure 
relief valve so that an initial damping rate, is replaced by C u l q2 

abo'-e a set motion rate, Y’^. 

Given the slope rate data, Y , far blade n at the damper location, c, the 

n 

c >. j 

initial damping rate is defined by |V I < t|. The damping is 

A 3r =- ( W (\) © <33 ‘ 


where 
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Beyond the pressure relief opening point the damping contribution is 
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The generalized force is developed for pitch horn bending and dynamic 
torsion. For pitch horn bending. 





(337) 


where M_ is the total feather moment as derived in Section 5.10. For 
rn 

the uncoupled dynamic torsion option. 
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in BLE axes. Since the blade elements loads are derived in BLn ax'is, 
the transform 
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is needed. The spring constant can be interpreted as 


K = M 
8 6 
PKn PKn 



( 3**0) 


where the generalized mass is computed continuously and u is the 

B PH 

nautral frequency of the uncoupled torsion mode, a program input 
constant. 


5.6.5 Quasi-static blade torsion . - To improve the pitch horn bending blade 
feathering representation a quasi-static blade torsion distribution is intro- 
duced. Quasi-static torsion is computed from the structural stiffness, 

at each station and the torque MX C at the shear center. The torque is 

s.mmed from the tip to the blade station in question as shown in Section . 3. 1* . 
Tht increment of twist produced at a blade station J can be displayed as: 
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assuming a first-order lag represents the torsional dynamics. The time con- 
stant is chosen to be representative of the blade first torsional mode 

frequency. 

To obtain this result, the available computation elements require some 
further operations. First, REXOR II conducts blade integrations from root 
to tip, in BLn axes. To obtain tip to root values: 
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Note the suosaation is conducted from root to the station J in question. 

Thus the j represents a summation vhereas the i represents a blade station. 



Second, these data are used to form the required torque at the shear center. 
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Small angles are assumed. The moments /R \ etc., act along the 

\ ^/BLn 

Bin axes and hence the matrix of lengths /X„„\ etc., are employed to 

V bt 7BLn 

obtain moments at the shear center which are then transformed into shear 
center axes, subscripted SC, parallel tc blade element center-of-gravity 
axes, subscripted BLE. 


The blade deflections and slope in BLn are also needed for the above 
expressions. 
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5.6.6 Quasi-static pitch horn bending. - To facilitate troubleshooting numer- 
ical instability problems as optional quasi-static pitch horn bending degree 
of freedom is available. The computation elements are the same as developed 
in Section 5.6.U except that the solution does not use generalized masses, is 
therefore an uncoupled mode, and is calculated externally to the main computa- 
tion flow. The formulation used is: 
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PH ’FnPH 


p FnPH 


PH 


(31*6) 


The dynamics are assumed represented by a first-order lag with B PH as the 

time constant. The variable is used to distinguish this formulation 

rnrn 

from the usual B pHn symbology. 

5.7 Shaft Axes Equations 


5.7.1 Transmission isolation mount. - The shaft equations couple the spring 
mounted transmissi' — , svashplate and rotor to the ground side of the mounting 
springs (fuselage). The fuselage is the reference coordinate set hence deriva- 
tives with respect to the rotor, hub, svashplate and transmission masses 
exist for the shaft axes equatic.s. 


5.7.2 Partial derivatives. - ry using rotor pse . oordinate masses only a 
few opera* Ions are required to assemble the couj. .ass terms using one 

derivative vector, |3 T q / 9 t s |- 
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This expression is a compact notation form of the development of section 5.5.3. 
Note that the angle to angle portion of the array is not a full transformation, 
but rather reflects the relation of a dependent coordinate to an independent 
Euler angle. 
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Hie hub to rv'ashplate partial i3t- /3t„1 is developed in the next section. 

I °SP H f 

5.7.3 Generalized masses . - The shaft axes matrix elements couple to the 
blade {A , B D „ J and swashplate generalized coordinates as well as to itself. 

Iflfl r HH 

Use is made of matrix notation and the rotor pseudo coordinate to produce a 
compact notation. 
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• Fuselage 


• Engine 


The six rigid degrees of freedom: X t I, Z, ♦, 8, ♦ are taken with respect to 

the stationary fuselage axes which are also the principal axes. The other 
elements considered are then referenced to the fuselage axes. The hub is sub- 
ject to shaft bending motions relative to the principal axes. The tail rotor 
is installed on the fuselage and rotates at the main rotor speed times an 
appropriate gear ratio. Positive rotations are defined as: 


• Hub 

• Swashplate 


same as main rotor 


• Tail rotor - Clockwise looking right 


• Engine - Counterclockwise looking forward 


The engine is treated as a rigid rotating body but the tail rotor is allowed 
to flap (teetering hinge, etc.). This flapping is considered secondary and 
enters only into the aerodynamic computations. The main rotor is allowed 
feathering, bending and twisting. 

5.8.2 Partial derivatives . - Elements used the reference axis masses and 
forces can by in large be conveniently related to either the fuselage or hub 
coordinates. Since the reference set is taken to be the fuselage coordinate 
set, no partials are required in this instance. 

Partials ’•elating hub coordinate motion to reference generalized follow the 
scheme given in section 5.5.3. 


where 
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The svashplate system is physically connected with the huh structure, .and 
partial derivatives describing its motion are taken through this intermediate, 
hub, coordinate. Note due to a parallelogram linkage the svashplate vertical 
motion is assumed to be unaffected by tilt angle. 




(357) 


5.8.3 Generalized masses . - Use is made of matrix notation and the pseudo 
rotor coordinate to express the reference generalized masses given in Table 4. 
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The mass of the fuselage is considered to contain -he engine and tail rotor 
masses, although the moment of inertias is treated separately. 

5.8.1+ Generalized forces . - The leads Lssociated with the six reference axis 
degrees of freedom are listed in Table 5. The tail rotor and engine are 
assumed to have shafts parallel to the fuselage axes. The transfer of the 
zerodynamic loads from tail rotor axes with origin at hub center and parallel 
to the fuselage reference axes is shown in the table. The fuselage aerodynamic 
loads include tail rotor and propulsion terms. Further development of the 
main rotor blade component loads is in Section 5.6.1 and the aerodynamics for 
all rotc-s and f ix^ ' surfaces is left to Section 6. 
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TABLE 5- - Concluded 



The angular velocities and accelerations associated with the engine and tail 
rotor require special consideration. Here the terms consist of a reference 
motion plus the turning due to the geared main rotor rate. 

Using the rate and acceleration Euler transforms for zero Euler angles 
(Section U.4.3): 
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5.9 Swashplate Equations 
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5 . 9 .I Partial derivatives . - The swashplate partial derivatives are readily 
obtained from Section 4.3.5- Using matrix notation 
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Since the swashplate axis is directly referenced to the principal (hub) set, 
the above derivatives are complete. The lack of translation to angular de- 
rivatives is explained by the parallelogram linkages used with swashplates 
to isolate the collective and cyclic inputs. The terms left out of the 
matrix indicate that the swashplate does not have a yaw degree of freedom. 


The reader should be aware that the angular notation 4>, 9, and V have two 
meanings, depending on whether they are in the numerator or the denominator 
of the partial. The numerator is the displacement of the mass element with 
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respect to the hub axis, whereas the denominator is the degree of freedom 
incremental variable. 

Swashplate motions pick up large inertia loads from the rotor due to blade 
feathering. Partials relating feathering to swashplate motions are 
assembled by first relating the feathering motion in the rotating system 
with feathering in the stationary system: 
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From Section 4-5-8, equations relating swashplate motions to the stationary 
feather angles give 
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where 


Also 
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The I T.k I matrix does not follow the conventional Euler angle notation 

L phJ 

since a desire existed to define as the angle the pitch horn to 

rH 

pitch link attachment point leads the blade. The overall derivatives can 
be put together as : 
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These partials are the elements of the 




vector. 


5.9.2 Generalized masses. - Table 6 presents the generalized masses which 
couple the swashplate motions with one and another and with th>? blade fuselage 
degrees of freedom. The table uses summations of the blade that arc described 
in detail in Section 5.6.1. 
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TABLE 6. - SVASHPLATE GENERALIZED MASSES 




Note p, q terms are the same for R and NR systems. 
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The moments used in these formulas are developed below. 

The feathering moment, M^, is taken to be composed of blade and friction 
loads. 
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The detailing of feathering moments due to blade loads, is accom- 


plished in Section 5.6.U. The friction load, follows the func- 
tion shown in Figure 30. By reducing to near zero, stiction 

rQ|DA. 


is obtained. Otherwise, if '♦ „ is large, the ratio 

rxl 9 JdK. 

determines the amount of viscous friction. 
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The remaining portions of the generalized force are the potential energy 
and dissipation functions. First consider the angular potential energy 
terms which model the swashplate tilt spring rate. This spring rate has 
a center dead-band, an operating range spring rate, and a high spring rate 
to simulate a travel limit stop. 

Consider the normal operating range spring rate first. The swashplate 

springs are defined in control axes (Figure 31) as K, and K (can 

SP 0 SP 

be unequal in size). To find the elastic spring loads, the swashplate 
motions are first found in control axes as: 
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The geometric interpretation of ♦ is shown in Figure 31. 

Taking the swashplate deflections in the control axes, subtracting control 
inputs and 0 C> and using the inverse transform, the swashplate spring 

terms in swashplate axes become: 
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where the subscript (l) is used to distinguish these values (used in sub- 
sequent logic calculations) from the final expressions developed below. 

Substituting for the swashplate motions in terms of the swashplate axes 
and rearranging: 
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where 
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and 


w- 


‘K cos 2 * + K sin 2 * 1 (K - K ) sin*, cos* 

♦sp c 0 sp c ♦sp 6 SP C C 


( K * ' K 0 ) 
L\ sp sp/ 


+~ — 


2 2 

sin* cos* K. sin *_ + K. cos *„ 

C Cl Tpn v C 


r SP 


SP 


Note: [K sp ] is a symmetric matrix of constants. 

The center dead band is modeled by the following logic. 
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(392) 
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(396) 
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is added to 



SP ) 


(398) 


to account for a limit travel stop. The limit deflection for the 
swashplate is 
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where 6 is the circular stop swashplate deflection limit. 

b »b.r 

The angular damping term is analogous to the spring load: 


where 



' 3B ] 


r ■ 
■4 

V SP 

V SP 


3B 

■ * m- 


®SP 

39__ 


SP 




has the same formulation as 




( 1 * 00 ) 


Control friction is treated as having rotating and nonrotating components. 
The rotating component has already been discussed as part of the feathering 
moment. The nonrotating component is applied to the swashplate. It has 

the formulation shown in Figure 31 with a change in labels such that 


is either * gp or e $p . and is either or 


The vertical potential energy term is described as: 
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Otherwise 


* K \ + K 2Z (z - ^ )♦ F ( 

3Z SP 1z sp 1 sp 2Z SP sp X SP ' 


(U02) 


F is a constant to center the gyro springs. 

v 

The spring rate is taker, to he K_ out to deflection and K 

Z SP SP SP 

beyond . 

A simple coupling from the rotary dampers gives the vertical dissipation 
function. 


3B • • • 
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where R^, R ?Q are coupling ratios. Note the effect of vertical 

motion on the swashplate tilt loads through the rotary dampers is assumed 
zero. 

To correlate with flight test records and/or to force the swashplate 
vertical response to cross the spring rate changeover a force offset con- 
stant is used. Introducing this constant into the swashplate vertical 
degree of freedom equation line, causes the variables, primarily the 
swashplate vertical motion, to shift and rebalance the equation. 

5.9.1* Control inputs. - The swashplate input is controlled by the pi lot' r 
cyclic stick. The input torques are: 
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The inputs are aligned with the control axis (Figure 31). 

Note the equivalence of forms in terms of angular commands 4 >p, 8 , 0 , or 

C V 

longitudinal stick (aft) Xo and lateral stick (right) Y c . 
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The controls are frequently linked to the svashplate through actuators which 
as a first-order approximation, can be simulated by a first-order lag. See 
Section 7.2. 


5.10 Engine Equations 

5.10.1 Rotor azimuth and rotation rate . - The program allows a variation of 
rotor speed in maneuvers due to variations in the torque required by the 
various rotors and in the torque supplied by the engine. The dynamic system 
rotates as a rigid, geared unit. That is, the shafts are not allowed elastic 
windup. The main rotor speed, ^ and hence the engine speed, is referenced to 

the fuselage, and not to inertial space. The displacement is the azimuth 
of the number one blade. n 

5.10.2 Engine model . - Figure 32 illustrates the engine model used in the 
program. The figure also plots typical engine torque characteristics. The 
model represents the first-order lag power response characteristics of the 
free turbine powerplants commonly used in rotorcraft applications. 

Being a perturbation model, the engine is referenced to its trim position. 

The change in engine torque in a maneuver is 


M. 
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( 1 * 05 ) 


where 0 S My. i M y . . The zero limit occurs if the overrunning 

ENG ENG, MAX 

clutch disconnects the engine in the transition to autorotation. The maximum 
value corresponds to the engine shaft torque limit. 


The gas generator, speed, is a degree of freedom. It is considered 

a secondary degree of freedom in that the coupling through the generalized 
masses with the primary degrees of freedom can be neglected. An equation for 
the generation speed can be supplied from its torque characteristics: 


I GEN ^GEN + C GEN V GEN = “ K ENG1 V ENG * K ENG2rENG ~ V ENG, TRIM , 


(406) 
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The terms on the left represent acceleration inertia torque and steady- 
state torque. On the right, the fuel control causes torque to he added if 

the engine speed drops below the trim value. The term exists since 

the control is modeled with simple lag. Restating this equation, using 
rotor speed and a generator time constant, gives: 
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where %> 
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is the order of a second. 


The engine droop characteristic can be used to size the engine constants. 

• • mm 

With = ♦_, = 0, substituting the generator equation .ne engine 
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equation and rearranging. 
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Only A incremental changes are of interest. The bracket subscripted R 
indicates the torque is determined at the rotor speed and includes the 
engine gear ratio. The term on the right is the static droop line shown 
in Figure 32. This plot also geometrically interprets the partial 
derivatives on the left. 

The generator speed ^ is not given a reference. Its value is zero 
when trim is completed . 

5-10.3 Partial derivatives . - Shaft rotation not only involves blade root 
rotation 4^, but also feathering notions. The feathering partial is 

obtained by differentiating the feather angle equation in Section U.5.S: 
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A set of partial* are defined to relate the various rotating components to 
t*!» rotor shaft and to the reference set. 
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(Ml) 


( 1 * 1 2 ) 
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5.10.1 Generalized masses . - The engine degree of freedom couples with 
every other degree of freedom. Equations for the engine generalised masses 
are given in Table 7. Matrix notation is again used for compactness. Note 
the transmission is modeled as a non rotating mass, and therefore dees not 
appear in these masses. The engine degree of freedom contains not only 
rigid body motion of the rotor blades, but also blade feathering. The 
feathering contribution is a minor contributor for some of the mass matrices 
and has been neglected. 
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5-10.5 Generalized forces . - Only one generalized force is needed. 


F n (Fiii)jijp " GgjjQ 1^ (Fp ” ^pur. 




F ENG R 


°tr (F * TO ' r iy ( *F + g tr V> 

A 


- (I ♦ I ) ^ + 
^ zz sp R 


Nb 


n=l 




( 421 ) 

The main rotor contribution (inertial and aerodynamic loads) is eiven in 
Section 5-6.U. The tail rotor aerodynamics are described in Section 6.5. 
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6. AERODYNAMICS 


6.1 Introduction 

Other than gravity, the external loadings acting on the REXOR II equations of 
motion can be traced to aerodynamic sources. The following subsections trace 
the source, nature and use of these aerodynamic loads. 

6.1.1 Aerodynamic forces producing surfaces considered . - The aerodynamic 
loads considered in REXOR II are divided into the categories of (l) associated 
with the main rotor, or (2) the rest of the rotorcraft (nonrotating surfaces 
and tail rotor). In view of the stated objectives of REXOR II, the program 
development emphasis is on the main rotor which is considered in Section 6.2. 

The nonrotating components consist of the fuselage, wing, vertical tail, lower 
horizontal tail, upper horizontal tail, tail rotor, auxiliary thrustors, mov- 
able surfaces on the wings and empennage, and dive brakes. Wake effects from 
the main rotor and wing are addressed in Section 6.3. The nonrotating load 
elements are mostly developed and assembled in Section 6.1*. The tail rotor 
equations, in integrated form, are developed in Section 6.5, and the auxiliary 
thrustor formulation is in Section 6.6. 

6.1.2 Use of forces generated . - As mentioned, the aerodynamic loads are in 
essence the external forcing functions of the equations of motion. Generally 
the developed loads are in the axis of the apparent air velocity of the loaded 
element. Thus transformations are required to put the loads into the refer- 
ence axes of the equation of motion considered. 


6.2 Main Fotor 

6.2.1 Overview . - To generate a main rotor model with sufficient detail to 
do dynamic investigations, a reasonably good quality aerodynamics presenta- 
tion is required. Tb this end a table lookup of blade section properties, 
multifunction inflow model, quasi-steady aerodynamics, and dynamic stall are 
used in REXOR II. 

6. 2. 1.1 Blade flow field . - As developed in the following subsections, the 
instantaneous blade airflow is the inertial velocity of the blade element. 
This velocity includes the motion of the principal reference set and the 
motion of the blade element with respect to the principal reference set. The 
calculation assumes the airmass is at rest, which is reasonable for dynamics 
investigations . 

6. 2.1. 2 Air pressure and angle of attack . - The dynamic pressure used for 
these calculations is based on sea level standard density. The loads are 
ratioed to the actual air density. 
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The angle of attack is the sun of geometric pitch angle and the instantaneous 
air velocity. The rate of angle of attack is also calculated and used for 
the transient blade aero loads. Sections 6. 2. 3. 3 and 6.2.3. 1 *. 

6. 2. 1.3 Forces and moments produced . - The steady blade loads are produced 
fron the air velocity components of Section 6.2. 3.1 and the coefficient *ata 

(C L 

developed in Sections 6.2. 3.3 (quasi-steady aerodynamics ) and 6.2.3. 1 * (dynamic 
stall). 

6.2.2 Concept of rotor inflow model . - The main rotor inflow model used ir. 
REXOR II is based on the air flow incident upon the rotor disc plus the air 
velocity imparted due to momentum exchange due to integrated blade span load- 
ing. This is to be contrasted with a formulation which tracks the rotor 
blade positions and the attendant trailing vorticies. 

The incident air flow is the inertial velocity of the rotor coordinates, and 
is directly available from the preceding mechanical development. However a 
number of assumptions need tc be stated and utilized to arrive at the induced 
velocity component of the inflow model. 



1. Only the vertical downvash and its variations radially and azimuth- 
ally over the rotor disk are considered. Induced swirl and lateral 
downvash components are neglected. 

2. Downvash effects due to unsteady aerodynamics are not treated here 
as an overall effect, but as a blade segnent condition in Sec- 
tion 6. 2. 3. 3- 

3. Rotor-induced flow distribution in hover and forward flight is pat- 
terned after Reference 6. This reference assumes a uniform loading 
in hover. Figure 33, from Reference 7, shows this distribution 
compared with typical loading and a triangular loading model. Fig- 
ure 3*» from Reference 6 shows the theoretical induced velocity dis- 
tribution in forward flight as a consequence of a uniform hover 
distribution. Thi3 data is fitted to slopes or a longitudinal skew 
as a function of speed in REXOR II. Lateral distribution remains 
uniform in accord with Reference 8, which corrects the lateral dis- 
tribution work of Reference 6. 

U. A variation in lateral and longitudinal induced velocity is included 
to account for roll and pitch aerodynamic shaft moments. 

5. Lifting line theory correction is accounted for by an effective 
rotor radius, BR. 


, Cjj, C^) of Section 6.2.U. The transient lift and moment effects are 
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Figure 33. - Blade loading distributions in hover. 
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6 


Root cut out effects are ignored 


7. Transient effects are simulated by a single time lag. 

6. 2. 2. 2 Steady state values. - The starting point for determining the down- 
wash is momentum theory as applied to an elementary dA: 

dT * (#) - (% m “) *1 

= p \/ ♦ V H 2 + (v H - v t ) 2 dA 2 v^ (J* 22) 


The thrust increment la dT, dm/dt is the flow of air through the rotor 
disk with resultant velocity p is the air density and the down- 

wash velocity. The velocities are taken in hub coordinates and no effort 
is made to account for rotor tilt. 

The thrust expression above is used to define the following induced 
velocity components. 

Average component, w^ 

Longitudinal variation with pitching aerodynamic 
moment, q iMR 

Lateral variation with roll aerodynamic moment, p.. ro 

lMn 

The downwash velocity becomes 


W i S W iMR + r ^iMR C0S + r p iMR Si ” *R (U23) 


The coefficients can be evaluated by equating the thrust and moment values 
for the main rotor equations to the integrals of the momentum expressions 
at hand. First consider the thrust expression. The evaluating task can 
be reduced by employing some boundary conditions. For rotor thrust only 
(no moment), q^^ and pj™ =0. A convenient expression for the elemen- 
tary area, dA, is shown in Figure 35. While radial anulii would serve 
for thrust integration, the form selected is particularly suited for the 
moment expressions. 
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For the average rotor thrust 



( 1 * 21 *) 

A further assumption is required to solve the square root of this expression 
and the corollary momentum equations. 


For forward flight 


w. « 
1 



w. 

iMR 


) = (constant) 


( 1 * 25 ) 


Completing the integration gives 



P7r(BR) V iMR 2v 


iMR 


(1*26) 


Next consider the case of no rolling moment; i.e., only thrust and pitching 
moment. Figure 35 is used with the incremental strip considered to be right- 
left oriented so that all equal values of q are integrated at once. 

Then, 



(1*27) 
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or 



r ■IBB )' 1 
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iMR 


q iMR 


(U28) 


Likewise for rolling moment, and using fore-aft increment strips gives 




2 P 


iMR 


(U29) 


Note the subscript A on F , , and VL denotes the aerodynamic 

mR,H ~MR,H MR ,H 

component only of main rotor loads in hub axes. 

The foregoing expressions are now developed for hovering and low-speed flight. 
In this condition. 



+ v„ 


« 



(1*30) 


Integrating gives 


H =P ’ (BR)2 
Ma,n 
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iMR 


(«31) 


and 
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The consequence of cyclic, first-harmonic down wash was explored in 
Reference 9. Their conclusion, which parallels Lockheed's experience, is 
that the phase and magnitude of the flap response of a hingeless blade to 
cyclic feathering is markedly affected by cyclic downwash. The shaft moments 
variation with feathering angle and the phase angle between flap and feather- 
ing are both reduced with cyclic downwash, the effect being greater in hover 
than in forward flight. 

A physical interpretation can be rationalized for the formula above, at least 
in hover, in that the aerodynamic thrust and moment produces a flow of linear 
and angular momentum. Imagine the flow as a continuous stack of disks having 

2 

"iass per unit thickness P*(BR) and dimetral inertia per unit thickness 
o (* (BRj^/li ) . 2 Wjjjp* 2 p iJfR and 2 q 1MR are the final, far downstream posi- 

tion, values of induced velocities obtain by these disks oriented with the 

flow. The terms p«(BR)^V. anv. d*(BR;**A V. are the mass flow per unit time, 

xi 

and the moment of inertia flow per unit time through the actuator disk, 
which times 2v^ MR , 2p^^ or 2q iMR is the gain of momentum. 

For programming purposes, an empirical blend of the forward flight and hover- 
ing sets of expressions is used. The limiting cases of the empirical set give 
the derived cases. The expressions used are: 
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(1*35) 


(1*36) 


6.2.2. 3 Variations in forward flight and ir g round effects . - The previous 
development can be assembled and combined with linearized forward flight dis- 
tribution and ground effect factors. 

= W iMR f iMR f 1 + K iMR E cos( *R + ^BLn + V] 


+ rp iMR Sin ^R 


W + rq iMR C ° S( *R + ^BLn } 


(1*37) 
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The ground effect factor, f .^, and the longitudinal linear gradient factor, 
^iMR* accounts ^ or forward flight. 

Note this formula is in rotating coordinates, and that the forward flight dis- 
tribution actually is applied along the line of the apparent airflow, In 

doing this the distribution is valid for forward flight, sideward flight and 
sideslip conditions. The angle 1> v is 


r 


w 



(1*38) 


as shown in Figure 1*5. 


The aerodynamic moment factors, 
axis. 


and p.^, remain attached to the hub 


The downwash factor K.^ , as explained 
of the wake angle defined as 


in assumption 3, 


is given as a function 


X iMR 


Tan 
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/ 2 2 

^ * V H 

W iMR ~ W H 


(1*39) 


which is ze’-o in hover and rear 90 degrees in high-speed flight. The function 
can be constrained by a number of factors. In hover, the value is zero. A 
90-degree value of about 1.6 can be read from Figure 3 1 * - Also from this fig- 
ure a set of linearized distributions is read, and plotted as Figure 36. 

The ground effect factor. 


f = 1 _JL/_EA i_ . 

iMR x 16 \ h / 2 2 

, . “H + V H 
1 T — 

w iMR c 


( 1 * 1 * 0 ) 


is taken from Reference 10, where h = -(Z0„)_ is from Section U.5.1. 

n £ 

6. 2. 2.1* Downwash transients . - Downwash transients exist due to an apparent 
mass associated with the induced flow field. Work by Peters, et al. derive 
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the expressions for collective and cyclic dovnvash including unsteady 
components. A good suaaary of this vork is Reference 11. 

Converting the referenced vork to dimensional form gives equations comparable 
to (U3U, U35, U36) . 
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These differential equations are solved for v.._, P m and q.,_ using numeri- 

. i „ , \ . » lMn iMn lMn 

cal (Euler) integration. 

6.2.2. 5 Iteration of dovnvash solution . - As is the case with any rotary-wing 
loading calculation, there is an interplay between the dovnwash variance from 
calculating the loading and a variance in the loading from recomputing the 
dovnvash. A common practice is to solve an iterative loop to satisfy both 
equations (i.e. , lift and momentum). In REXOR II the iteration does not take 
place independently, but proceeds stepvise with the rotor azimuthal advance. 
With the normal, rapid convergence of the iteration the solution will essen- 
tially be complete with the step advance. Hovever, large step sizes vill 
incur an additional dovnvash time lag. 

6.2.3 Blade element velocity components . - In the following subsections the 
blade aerodynamic loading is categorized and developed along two lines. They 

are: 
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• Steady-state aerodynamics 

a Transient phenomena consisting of quasi-steady aerodynamics and dy- 
namic stall. 

6. 2. 3.1 Sources and resolution from Made motion . - The steady aerodynamics 
are based an the air velocities while the quasi-steady aerodynamics (from 
flutter theory) and dynamic stall depend on accelerations . 

The air velocity is the hlade mechanical velocities suaaaed with a component 
due to dovnvash. In a similar manner, the air acceleration is taken to be 
the mechanical blade accelerations minus the dovnvash accelerations. The 
dovnvash formulation as developed in Section 6.2.2 allows for lags, and it is 
these lags that result in dovnvash acceleration terms. 

6.2. 3. 2 Steady aerodynamics . - The air velocities relative to a blade section 

are desired for an axis system with origin at the quarter chord to match the 
airfoil table data. Frcm Section 4. 5- 5. the mechanical blade velocities rel- 
ative to the free stream or earth axes are available as { Y BLE* ^BLeP* 

The desired relative air velocities at the quarter chord (or blade Bln refer- 
ence axis) are . 
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Where the second vector on the right is the dovnvash velocity developed in 
Section 6.2.2, and the third transfers the velocity from the BLE reference 
point at the blade center of gravity back to the quarter chord. The distance 
i3 positive with the center of gravity ahead of the quarter chord. For 

rotational convenience. 
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The angle of attack is defined as 


°i/i» . SIn ' 1 >*•«) 

Airflow aspects of quasi-steady aerodynamic formulation are developed at this 
point for convenience. The quasi-steady aerodynamic contribution is conceived 
as composed of circulatory and noncirculatory components. The circulatory 
components are taken to be equivalent to finding the aerodynamic force, and 
moment coefficients are based on an angle of attack at the three-quarter 
chord: 
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3/ * c 
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As such, the effect of angular rates is included in deriving the steady aero- 
dynamic coefficients. The formulation above does not attempt to account for 
local dovnvash rotation or curvature and its chordwise variation. The net 
result is that aerodynamic coefficients determined in Section 6.2.1* are com- 
puted with c - 

A number of quantities used in the dynamic stall computations. Section 6.2.3. 1*, 
are also available from the previous mechanical development . They are also 
defined here for convenience. First, the angle of sideslip appears only in 
the dynamic stall formulation. For this purpose it is defined as 


A = Tan 
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(1*1*8) 


Also, dynamic stall is based on the time derivative of the angle of attack at 
the three-quarter chord: 
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The inclusion of the gravity tern places these accelerations in a true iner- 
tial axis systea, not earth inertial axes, as appropriate for aerodynamic 
calculations. See Section U. 5-1. Gravity does cause buoyancy forces, but 
these can be ignored. The turning acceleration components are also subtracted 
• to produce linear accelerations which correctly aodel the blade element inci- 
dent airflow. The gravity vector can be obtained from hub values as: 
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By differentiating the down wash velocities, the down wash acceleration is 
obtained: 


(“ BLE W.Btn) 


. . ( VaLE ) W. 
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6.2. 3.3 Quasi-steady aerodynamics . - Quasi-steady aerodynamics is accounted 
for in REXOR II by incorporating the terms frcm the two-dimensional flutter 
theory of Theodorsen (reference 12). In the REXOR II analysis, Theodorsen's 
lift deficiency function C(k) is taken as unity. This means that the flutter 
theory presently incorporated neglects shed wake effects, or in physical terns 
does not account for the phase change between blade element lift (or pitching 
moment) and angle of attack, due to shed vorticity, or the assumption of 
quasi-steady aerodynamics is expressed by C(k) = 1. 

Referring to a classic text on aeroelasticity by Bisplinghoff , Ashley, and 
Hofftnan (Reference 13), the expressions for lift and pitching moment are 
given as: 
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In REXOR II t the blade aerodynamics and quasi-steady aerodynamics are 
referenced to the local section quarter-chord properties . This is done 
because the majority of available airfoil data uses this reference. Rote 
that the final aerodynamic loads are translated to the local BUB axis (c.g. 
location) for use in the equations of motion. 

Reviewing the above expressions, and referencing the rotation point to the 
quarter chord gives a * -1/2. If we take C(k) as unity, replace 2 « for cir 
culatory lift by (dC^/da), and substitute c/2 for the semichord b, these 
• equations becose 
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and: 


M-if- 
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Rote that the entire last term in the moment equation vanishes with a = -1/2. 
Referring to the lift expression, noncirculatory aerodynamic lift is accounted 
for in REXOR II by the first term in which h ♦ Ua are combined into il^in 

blade elment coordinates. The second term results from table lookup where 
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in which the angle of attack is previously computed from 


“3/Uc = 



(1*58) 


The a within the brackets is identified as 9 , the actual physical angle of the 
blade with respect to the freestream direction. The a on the left hand is 
that due to the air velocities which include the plunging velocity h and ro- 
tation component c/2a. Hence p T * a and p_ T _ * a. 


BLE 


The total aerodynamic pitching moment is the sum of the quasi-steady loads 
computed above and the table lookup blade section properties (Section 6.2.4). 
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6. 2. 3. 4 Dynamic stall . - Dynamic stall is included in REXOR II based upon 
the Boeing-Vertol formulation set forth in References lU, 15, and 16. It is 
similar to the treatment of dynamic stall in the Bell C-8l program. A com- 
parison of REXOR II with the C-8l program is given in Appendix IV, pages 393- 
U04, of Reference 3. Dynamic stall is specifically addressed with respect to. 
the two programs beginning on page 395 of that report. A significant point 
of difference between the treatment of dynamic stall in the two programs is 
that C-8l puts a 20-percent limit on the angle-of-attack overshoot in obtain- 
ing the dynamic maximum lift coefficient, whereas REXOR II has no limit. The 
correctness of the treatment of dynamic stall in either program is difficult 
to assess since the concensus of researchers in this area is that current 
methods are empirical at best, and much research still remains to be done in 
this area. 

Reference 14 notes that, "The trends show that compressibility effects reduce 
dynamic-stall delay, and at about M = 0.6 no dynamic-stall delay is evident." 
For this reason an upper Mach number limit of 0.6 was implemented in the dy- 
namic stall calculations for REXOR II. The test data obtained by Boeing 
Vertol and given in the references cited was for the Mach number range 0.2 to 
0.6. As implemented in REXOR II if M < 0.25, the value M = 0.25 is used in 
the analytic expression for developing the stall hysteresis loop. 

Reference 15 notes that it was found that, "airfoils used currently by the 
helicopter industry had stalling dominated by leading edge stall. For this 
type of stalling process, the dynamic C extension was proportional to the 

Xj 

time rate of change of the angle of attack." 

In that reference, so as to use static airfoil data as much as possible, 
static stall and dynamic stall are empirically related by developing a refer- 
ence angle of attack given by 


°REF 


a 


Sign( “ } ) 


(459) 


in which. 


Y = log e 


0.601 

M 


(460) 


and is physically related to dynamic stall delay. a is identified as 
and 3 ' kc 


M = 



(461) 
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As noted in Reference 16 in regard to dynamic stall ..." as a blade element 
reaches and exceeds the static angle of attach, stall does not occur as long 
as a sufficient, positive time rate of change of the airfoil angle of attack, 
o, is present." The experimentally derived equation for dynamic stall delay 
is given in the reference as 


dynamic stall delay * y 



(462) 


where 


ca 

2V 


k 


(U63) 


the blade element reduced frequency. 

Referring to the gamma expression, we note that T + 0 as M > 0.601, which is 
the upper limit for Mach number values for dynamic stall calculations. Also, 
note that Y 1 as M -*■ 0.2211, which is approximated by the value of M = 0.25, 
the lover limit in REXOR II for dynamic stall simulation. 

The termOj^y given above is also called the dynamic angle of attack (Refer- 
ence 15) and given by the notation <* . 

6.2.3-l*.l Lift accounting for dynamic stall . - Using the reference or dy- 
namic angle of attack computed from o^, the REXOR II program implements the 

"Fast Aerodynamic Table", Section 6.2.1*, subroutine and determines the lift 
coefficient, C^, corresponding to a and the freestream Mach number for the 

specified blade element and blade azimuth position. Also computed at the 

given Mach number are the C T for zero angle of attack and the C T for a small 

L L 

increment A a with respect to zero. Yawed or radial flow is accounted for by 
computing the yaw angle of the flow given byr 


A = 


Tan 



(U6U) 


where U c and U n represent blade spanwise and chordwise components of flow 
respectively. 
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The slope of the lift curve is then found from: 



C L (a REF. M) - C L (0 » M) 
“REF cos A 


(U65) 


It can be argued from physical reasonings that the dynamic lift-curve slope 
cannot exceed the static life-curve slope. As a check, REXOR II also 
calculates : 



CjAa, M) - C l ( 0,M) 
Act 


( 1 * 66 ) 


Only in the event * s 8 reater than (aC^/Sa)© ^ is the latter value 

used to calculate C . Otherwise C is calculated by 

L L 


a+C L (0,M) (1*67) 



The ability of this approximation to describe mathematically the lift hyste- 
resis characterized by dynamic stall is shown in Figure 37, which compares 
analytical results with experimental two-dimensional airfoil data. (From 
Reference 16. ) 

The component of the lift force per unit span acting normal to the blade 
chord axis and including dynamic stall effects is then calculated from 


PV 2 

af no = C L c T coso 


(460 ) 


The total normal force is determined by adding to this term the drag compo- 
2 

nent, c PV /2 sin a , and the unsteady aerodynamic terms discussed in the 

previous section. To account for dynamic stall effects on drag, two- 

dimensional drag coefficient data are used, but as determined at a , not a . 

REF 

This is consistent with Reference 15. 
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The component of the lift force per unit span parallel to the blade chord axis 
is found correspondingly from: 


C- c PV 2 

AF C = — 2 3ina 


( 1 * 69 ) 


The total chordwise force is then obtained by auding the corresponding drag 
coefficient term multiplied by cos a . 

6. 2. 3-1*. 2 Pitching moment accounting for dynamic stall . - For determining 
pitching moments due to dynamic stall (see Referer ce 13), the reference or 
dynamic engle of attack given °y a RE p must be modified. In REXOR II, this is 

accomplished by multiplying the second term by an empirical constant, K. 
Hence, 


a REF = °DYN a ~ K ( | 2V j si 8 n ( a ' j i^lO) 

K is selected based upon the dynamic stall characteristics of the airfoil. 

In general it has been found for conventional rotor blade airfoils that K 
should be selected so that 


a » 
REF 


a 

REF 


+ A a 


(U71; 


where da is of the order of 2.5 degrees. With a’ calculated from the above 

nr*r 

equation, the moment coefficient is determined from tables such that. 


C M = C m"W H) 


( 1 - 72 ) 


A comparison of test and theoretical dynamic C from Reference 15 is shown in 
Figure 36. 

The total pitching moment acting per unit span on a blade element is then 
given by: 


T(i) = 


C M c2£ T~ ~ F N0 S Y U) + 


~quas 4 - steady’ 
aero 
terms 


(Section 6.2. 3. 3) 


( 1 * 73 ) 



ANGLE OF ATTACK 


Figure 38. - Dynamic stall - moment coefficient vs 
angle-of-attack hysteresis loop. 




vhere S^(x) represents the distance from the aerodynamic center to the blade 

elastic axis, and the quasi -steady aerodynamic terns are included as described 
in Section 6. 2. 3. 3. 

6.2. U Coefficient table lookup - overview. - In cataloging blade section 
aerodynamic data, C^, and C^, there are two procedures available. 

• Curve fit *he aerodynamic data to the specific airfoil geometry being 
investigated for the range of Mach number and angle of attack to be 
considered. 

• Tabulate the data as a function of performance and geometric param- 
eters, and interpolate to the exact conditions at hand. 

REXOR II uses the second procedure. The data consists of C^, and tables. 

Each table is tabulated as a function of angle of attack and Mach number. The 
table format is organized identically to the Army C-8i program. Thus C-8l 
airfoil decks may be directly used in REXOR II. 

A table set of KACA 0012 section characteristics is included as part of 
RBCOR II. Two external tables may be used; the first of which overrides the 
resident 0012 data. Changeover of external tables occurs at a preselected 
blade radial station. 

6.2. U.1 Inputs and outputs. - Each table (C^, C^, C^) has a separate angle 

of attack entry and a common Mach entry. The separate entries are used for 

dynamic stall calculations. The outputs in addition to C , C and C u are the 

L D M 

zero angle-of-attack C and C T vs angle of attack slope. 

L L 

6.2.5 Blade element and rotor aerodynamic loads summary . - The required 
loads for use in the equations of motion are in BLn axis. Development to 
this form from BIE axis about this quarter chord point is covered in Sec- 
tion 5 .6.U. The BLE axis form is: 
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6.3 Interference Terms 

6.3.1 Nature of the Phenomenon. - In the process of producing lift, the var- 
ious parts of the rotorcraft impart a net momentum change to the air mass 
opposite to the direction of the force produced. This induced air velocity 
from the momentum change impinges upon other elements of the rotorcraft chang- 
ing their aerodynamic behavior. 
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The sources of interest are the main rotor and wing (or lifting body charac- 
teristics of the fuselage). The surfaces being affected are the wing plus 
fuselage and the empennage. The impinging velocity is expressed in Zp 
(fuselage vertical) axis as a percentage of the source flow and a function 
of the wake angle of this flow. 

A second interference velocity source is to consider the circulation part 
cf the Theodorsen function. Here the wing or wing equivalent of the fuse- 
lage is producing lift at the quarter-chord point according to the air 
velocity at the 3 A- chord location. Accordingly, the vertical component 
of air velocity at the wing includes a ccaroonent. 


i C wing q F (1*76) 


Here the wing quarter chord is assumed to lie on the Y-s- axis. This com- 
ponent is also effective at the horizontal tail via the wing to horizontal 
tail downwas'n factor. 

6.3.2 Rotor to wing/ fuselage. - The down wash function (percentage of source 
flow) used in REXCH II is a lookup table of dovnvash factor, F ( X ) and 

idealized main rotor wake angle 

where. 


*MR = tan 


\ H iMR / 


(1*77) 


rhe table data is linearly interpolated to the required wake angle value. 

The fuselage reference downwash velocity at the wing (or equivalent) then is 


! s V _ V W y 1 * — r n 

WING F iMR ^ ; MR-V 2 WING q F 


(1*78) 


and taking tine derivatives. 


i = v _ v v(\\ + _L r n 

WING F iMR V * ; MR-W 2 WING q F 


(1*79) 
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The total air velocity to the ving/fuselage is 


WING 


b 


1/2 


♦ v. 


+ w. 


WING 


and the angle of attack is 


3 WING 


tan 


-1 


(^) 


( 180 ) 


( 181 ) 


The total velocity in the fuselage XZ plane is used in the empennage 
computations : 




1/2 


+ v 


WING 


(182) 


6.3.3 Rotor to horizontal tail . - A dovnvash factor F(X) mr _ H t between the 

main rotor and horizontal tail is computed in the same manner as F(X)^_ W 

from the main rotor wake angle X^. This data in conjunction with the wing 

to horizontal tail dovnvash factor is used to compute incremental air veloci- 
ties at the horizontal tail. 

Evaluating the main rotor increment. 


w iMRj W iMR ( F(X) MR-HT “ F(X) MR-w) 

A increment for an upper horizontal tail is likewise generated: 

W iMR = w iMR ( F(X) MR-HTU ' F(X) MR-w) 

J- vJ ' 


(183) 


(131) 


6.3.1 Data sources . - The theoretical dovnvash factor ranges from 0 at X = o 
and 180 degrees to 2 at X =90 in the fully contracted rotor vake. Several 
sources of measured data are available to construct a distribution for a given 
configuration. Reference 12 gives isolated rotor data for field distances and 
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wake angle ranges suitable for F(X)^ and Reference 16 gives a 

good data set for typical wing locations. 

6.3.3 Empennage velocity components . - REXOR II nodels the empennage assembly 
as either part of the fuselage-wing aerodynamic table (tail on) or as a sep- 
arate set of aerodynamic loads (tail off). In either case a set of perturba- 
tion velocities is used. 

The wing to horizontal tail downwash factor appears explicitly as a quasi- 
unsteady aerodynamic term. An airflow time delay from the wing to horizontal 
tail is computed as 


At 



(W5) 


Using the downwash factor 3 e/3 a the vertical airflow component at the 
horizontal tail is 


At w 


wins 


3e ^HT j j ie 

3o " rflNG .. 


(U86) 


In a like manner the delay in sidewash gives rise to the term lyipA r 
on the vertical tail. 


Z 


^ytT./ 33 


The vertical incremental horizontal tail velocity is then: 


3c . 


Aw HT _W iMR I + l HT (%• + 3a VeNg/ 7 ^) 
+ ( T EL 6 EL " t HT 6 i )u. 


HT ' HT 


(1*87) 


The terms and introduce equivalent velocities due to elevator and 
horizontal tail incidence deflections respectively. 

Similarly for the upper horizontal tail: 


Aw t _, = -w.._ 
HTU iMR 


+ * 


IU 


HTU 


( q F + L *VING/ V X z ) 


( 1 * 88 ) 
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Assembling the vertical tail lateral incremental velocity: 



Where T pyp is used to introduce an equivalent velocity due to rudder deflection. 

A horizontal and vertical tail longitudinal total velocities are simply 
developed from a wake deficiency factor 


u 


HT 


u_ n _ 
F HT 


( 1 * 90 ) 


U HTU = U F n HTU 
and 

U VT “ U F n VT 


( 1 * 91 ) 

( 1 * 92 ) 


Then the total vertical velocity at the horizontal tail is: 


v 

HT 


_v F(x) 
iMR ' 'MR-V 


_ u e 
F FT 


♦ Av, 


HT 


and 


W HTU 


W F -V iMR F ^MR-W 


_u e ♦ 
F HTU 


Av, 


HTU 


(»* 93 ) 


(W 


The induced flow field angles and are a function of wing angle of 

attack, flap deflection and wing incidence change. 
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HT 



+ 


3e . 

— sin a 
da 


w 


. ac HT . ^ 3e HT . 

a °FL FL ^i w % 


( 1 * 95 ) 


C HTU 


3e 


HTU 


HTU 

3a 


sin a 


3e 
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HTU 
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The total lateral velocity at the vertical tail is 


v.„ * v_ - u_ o,_ + AV._ 
VT F 7 VT VT 


(U9T) 


Where is a sidevash coefficient from the fuselage . 


0 VT 


30 


sin e w 


(U98) 


6.U Body Loads 

In this section the aerodynamic contributions from the fuselage, wing and em- 
pennage are developed. These components are addwd together vitb the tail 
rotor loads. Section 6.5, and auxiliary thrustors. Section 6.6. A transfor- 
mation to fuselage axes is made from vind axes. 

The fuselage, wing and empennage data is composed of STATIC , DERIV, ar.d CON- 
TROL elements. The STATIC data are the steady state load components as would 
be measured in a vind tunnel. These data may be tail on or off. The DERIV 
data give additional loads due to velocity component variations from trim for 
tail on STATIC data as veil as steady offsets (unequal ving twist, etc.). 

Tail off DERIV loads use full tail velocity components rather than variations 
to generate the empennage forces and moments. 

The CONTROL loads account for flap, dive brake, ving incidence, aileron de- 
flection inputs via the control system. The rudder, upper and lover horizon- 
tal tail incidence, lower horizontal tail ele* ator inputs affect the empennage 
air velocity components, and are developed in Section 6.3.5* 

6.L.1 Ncnrotating airframe airloads . - The required loads are computed in 
REXOR II as the sum of steady-state forces and moments plus loads arising from 
stability derivative type terms and control surface inputs. The steady-state 
data are formed in terms of overall C^, C 

empennage assembly. 


D , and for the fuselage, ving, and 
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The static body loads are: 


where: 


F 'I 
Xt, 


f-C Q. " 

d t h a 

B 


I 
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0 

B 



% 




► - — 4 

► 

M.. 

°o 


a b 


0 
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M i q a c wing 

d-t— 


0 

5 

. J 

W, STATIC 

k - 


Q = — 

a 2 


PS, 


WIHG 



(1*99) 


( 500 ) 


The wing area, , and chord, , are actual or the equivalent of the 

lifting fuselage. 1 ‘Alternately, they nay be the reference length and area 
used for the available wind tunnel data. C^,, Cl_ and are linearly in- 
terpolated fron input data tables of C^, C^, C^, versus angle of attack, a . 
The data are interpolated on a„ T ,,_ from Section 6.3.2. The loads developed 
are in wind axis. 

The stability derivative lead contributions are computed as a 6 by 7 de- 
rivative matrix postmultiplied by a velocity component vector. For tail 
on fuselage aerodynamic data: 


23U 




(501) 


The components of the [F^j] matrix are discussed in Volume III. Two hori- 
zontal tails and one vertical are assumed. The upper tail is considered optional. 
The matrix also provides for asymmetric effects of wing incidence differential, 
for linear and quadratic sideslip variations and for wing roll damping. 


For tail off data the velocity vector is replaced by 



( 502 ) 
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These terms also produce forces and moments in vind axe 3. 

REXOR II includes the effects of flaps, ailerons and dive brake? in the non- 
rotating aerodynamic loads. The flap deflections are modeled as linear sta- 
bility derivatives of C^, C D and C^. The aileron load is the variation of 

aileron moment volume (rolling moment coefficient times wing area times ving 
span) vith aileron deflection. The input is for one aileron. Dive brakes are 
represented as a variation of drag area vith brake extension. The brake panels 
are assumed to be on the fuselage vertical axis and a distance - h^ below the 
fuselage reference. 

The desired loads are: 



The static and derivative terms are added to form the total body loads and 
transformed into fuselage axes. 
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-cos(B w ) -sin(B w ) 0- 

sin(B w ) cos(B w ) 0 

L 0 0 1 


and 


<=“(«„> ■ W 


sin<< V ' 'W /V x z 


cos(B ) = V /V 

X Z WING 


sin(S H ) = 


WING 


The air velocities v v , V , are defined in Section 6.J.2, 

WING X y T wtng 


6.U.2 Component additional airloads . - A total array, {qLOADs}, of 
rotor air loads is computed in fuselage axes. 
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The first component is described above. The tail lotor load vector auxiliary 
thrustor load vector are developed in the following sections. 


6.5 Tail Rotor 

A number of different levels of aerodynamic presentation accuracy and axis of 
representation may be used for tail rotor computations. In line with the stated 
objectives of REXOR II, a linear aerodynamic approach is used. A shaft axis 
reference is used for the analysis. In this system, the air velocity quantities 
involved are easy to visualise. Also, the flapping ard feathering motions are 
the true, measureable quantities. 

6.5.1 Formulations . - First, consider the airflow quantities available in 
fuselage axis. Figure 39. Note the tail rotor axes alignment with respect 
to the fuselage axis system. Formulating the components with respect to the 
fuselage axes gives: 


" u tr" 
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where 


a 

TR 


-rx 2 - Sin 8 
3 P w 


(513) 


The vertical component is approximated by the horizontal tail value. Terms 
subscripted w refer to velocity, w, and sideslip angle, g, at the wing. The 
flow effects induced by the wing-fuselage bombination are described by the 
vake velocity deficiency fact or n,^, the side wash angle, t tr , and its varia- 
tion with angle of sideslip. 

The velocity vector is then rotated through the tail rotor shaft lateral tilt. 
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Constructing the blade element tangential (U^,) and perpendicular (Up) components, 
as shewn in Figure Uo gives: 


u t ' (r o! tb * 5in *n> 

(515) 

U P * - t TB - V ITO - rS ' "tR 6cos »TB 

(516) 


where, 

is the tail rotor induced velocity 

and 


TR 


TR 


‘ W iTR 


(517) 


Expressing the blade element angle of attack as a small angle of 
approximation , 


a 


TR 



( 518 ) 


where. 



COS gl 


TR 


- sin* TR 
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A 6_ coupling is used to minimize tail rotor flapping, 
as a reduction in feathering for positive flapping gives 


Defining + 5 




0 • 0 ^^ 4* (5 a cos ^ 4- A b $ 

TR 7TR l na ¥ ?R 3TR °1_ ‘ *7R 

- R TR 

'Hie tail rotor analysis assumes no coning. 

The blade flapping, 6, is then 
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X TR 
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2l»0 


8 = - a, cos <i» 


(521) 






The tail rotor expressions of interest are the prime forces added to the 
fuselage system. First looking at the tail rotor thrust. For a blp.de ele- 
ment ve have 


d T, 


TR 


[■£ p a b c a U T 2 dr J 


(522) 


TR 


where a is the lift curve slope , b is the runber of tail rotor blades, and 
c is the blade chord (assumed constant). 


Substituting, 


d T TR = [*2 pa b c {6U T 2 + ^ U T )dr l 


J TR 


(523) 


Integrating for the entire rotor. 


T TR * h J »* » '/* / ( 8 U T 2 * U P °T> dr d j 


TR 


(52U) 


vhere B is the finite airfoil lift factor expressed as a so-called tip 
loss factor. 

Noting only even functions contribute to the integrand. 


T TR ~ 2 ° a(b C R) TR 


^®TR 3 ^ K) TR + 0 TR 2 “TR 


+ 2 (fiR) TR V TR + (ftR/ TR “TR *1^ T 6 3TR^ 


(525) 


Note the thrust is independent of the longitudinal flapping, but is a func- 
tion of lateral cyclic shown as lateral flapping times delta 3. 

Ihe required lateral flapping angle is obtained by equating the lateral 
flapping moment equal to zero. 
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gives 


' ' %R S3TB 


(527) 


To obtain the longitudinal flapping angle, the longitudinal rotor moment is 
formed and set equal to zero. 
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gives 
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In formulating the tail rotor drive torque, the blade profile drag is ex- 
pressed as 


C D = 


kC L 


(530) 


vhere C L is the average lift coefficient. Reviewing the thrust equation 
with a constant (average) lift coefficient gives 
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TR 


gives 


C L * 6 T TB /(‘"’t* A TF ( fl3 '“"’TO -I B "f 2 )) <532> 

The drive torque is expressed as the reaction to turning the tail rotor 
shaft. The plus sign is associated with a clockwise sense of rotation 
when facing a left-hand mounted tail rotor. 
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where 
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integrating 
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(continued on next page) 
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The remaining load tern in X^. Using the sane formulation methodology, 
Q ^ = - 


0 b c r D U T 2 sin* 
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Making small angle approximations , 

d = j-p(b c)^ [- C D U T 2 sin * 

+ a( 0M ? U p ♦ U p 2 ) sin* 
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integrating 
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The induced velocity is calculated from simplified momentum balance. 


W iTR * T TR P1,R B ^TR ^“TR + V TR + “TR ^ ) ^ 539 ^ 


Normally, the thrust, flapping, and induced velocity equations are solved as 
an iterative set. In REXOR II, these equations are solved for every pass 
(azimuth step) of the main rotor, and the tail rotor set convergence is assumed 
a priori. 

Note that the pitch- flap coupling does not appear in the expressions developed. 
This is due to the equivalence of flapping and feathering, coupled vith the 
absence of lateral flapping. 
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6.5.2 Airloads - control settings . - The force and moment terms are assembled 
for use in the overall fuselage loads. Section 6.fc The pilot control is the 
rudder pedals 0 . 



is a factor to account for in blockage on the tail rotor thrust. The 
equations for X^p, and Q^p are based on P q , the sea level density. The 
sign of G^p allows for a tail rotating in a negative direction; i.e., upper top 
moving forward. 

6.6 Auxiliary Thrustors 

REXOR II models the compound helicopter configuration by the inclusion of an 
auxiliary source of forward thrust. A perturbation bypass jet math model is 
used. It is assumed that all the thrust units are at the same setting. Further- 
more the dynamics of the engine rotating mass are ignored. 
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6.6.1 Formulations and airloads . - Based on a perturbation model the thrust 
for all units installed is: 
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where M is the freestream Mach number and <5 represents the total engine 
P P 

control parameter . 


Hie engines are located at height h p and distance aft of the fuselage axes. 
A thrust angle 0 q is also assumed. The engine contributions to the fuselage 
aerodynamic loads^are then: 
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7. CONTROL SYSTEM 


*» 


7.1 Overview 

REXOR II models vehicles ranging from pure helicopters to winged helicopters 
to compound helicopters with conventional airplane control surfaces. The 
control system is modeled as a set of pilot controls (stick, rudder pedals, 
collective, etc. ) which are coupled to the helicopter and airplane aerodynamic 
surfaces through a set of overall linkage factors (gains). These gains are 
slaved to a master control (phasing unit) which can he varied from the ex- 
tremes of pure helicopter to pu’-e airplane type of controls. 


7 . 2 Pilot Controls 


To simplify the operation of REXOR II the control inputs are mostly expressed 
as a percent of full scale (maximum input). The pilot inputs are: 
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c,P 


C,p 
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o,p 


R»P 

5 db,p 

S FL,o 

5 iv,p 


°iHT,p 


Longitudinal stick 

Lateral stick 

Rudder pedals 

Propulsion setting 

Collective blade angle 

Rotor speed setting 

Dive brake extension angle 

Flap extension angle 

Command wing incidence change 

Command horizontal tail incidence change. 


Pilot controls are combined with trim (T), initial condition (IC), rigging 
offset (subscript 0), and stability augmentation inputs (SAS). These com- 
bined inputs then operate the rotor and fixed aerodynamic surfaces. Scaling 
factors (K) convert the percentage inputs into angular and linear deflections. 
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( 5 * 47 ) 
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6 iHT = K iHTFS ( * 6 iHT,p + ^ 5 iHT,IC^ 

6 iHTU = K iHTFS ^iHTU.IC^ 
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The factors C- , G_ T , G , G , G.,.- are the slaved gains controlled by the 
C tL AJLL in nUiJ 

phasing unit. 

The quantities X* , Y* are processed through a first order lag and rate 
c c 

limiting prior to being applied as svashplate input commands, X and Y . 

c c 
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The fixed aerodynamic surface motions are shown in Figure 1+1. The pilot 
inputs are depicted in Figure 1+2. 


7.3 Stability Augmentation Systems 

REXOR II incorporates stability augmentation inputs to the lateral and longi- 
tudinal cyclic inputs , elevator, aileron, tail rotor collective, and rudder. 
These SAS inputs are all derived from fuselage axis angular rate information. 
Signal processing consists of a low frequency washout and limiter applied to 
all throughput. A first order lag is also used on some signals. 

The SAS coefficients are also computed on the basis of percentage of full 
scale deflection of the pilot control they are connected to. The same scal- 
ing conversion factors as used for the pilot inputs are applied to the SAS 
outputs . 

The six SAS channels are shown below. Figures 1+3 through 1+8. 
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Figure 4l. - Fx .. ju aerodynamic surfaces. 
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Figure 42. - Pilot controls. 
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